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^\ , Abstract. We analyse IRAS and COBE DIRBE data at wavelengths between 2.2 and 240 /im of the central 500 pc 

of the Galaxy and derive the large-scale distribution of stars and interstellar matter in the Nuclear Bulge. Models 
of the Galactic Disk and Bulge are developed in order to correctly decompose the total surface brightness maps 
of the inner Galaxy and to apply proper extinction corrections. The Nuclear Bulge appears as a distinct, massive 
disk-like complex of stars and molecular clouds which is, on a large scale, symmetric with respect to the Galactic 
Centre. It is distinguished from the Galactic Bulge by its flat disk-like morphology, very high density of stars and 
molecular gas, and ongoing star formation. The Nuclear Bulge consists of an R~ 2 Nuclear Stellar Cluster at the 
CIh' centre, a large Nuclear Stellar Disk with radius 230±20 pc and scale height 45±5pc, and the Nuclear Molecular 

Disk of same size. 

The total stellar mass and luminosity of the Nuclear Bulge are 1.4±0.6xl0 9 Mq and 2.5±lxl0 9 L©, respectively. 
About 70% of the luminosity is due to optical and UV radiation from young massive Main-Sequence stars which 
are most abundant in the Nuclear Stellar Cluster. For the first time, we derive a photometric mass distribution 
for the central 500 pc of the Galaxy which is fully consistent with the kinematic mass distribution. We find that 
the often cited R~ 2 distribution holds only for the central ~30pc and that at larger radii the mass distribution 
is dominated by the Nuclear Stellar Disk which has a flatter density profile. 
' The total interstellar hydrogen mass in the Nuclear Bulge is Mh=2±0.3x 10 7 Mq, distributed in a warm inner 

disk with J?=110±20pc and a massive, cold outer torus which contains more than 80% of this mass. Interstellar 
matter in the Nuclear Bulge is very clumpy with ~90% of the mass contained in dense and massive molecular 
clouds with a volume filling factor of only a few per cent. This extreme dumpiness, probably caused by the tidal 
stability limit in the gravitational potential of the Nuclear Bulge, enables the strong interstellar radiation field to 
penetrate the entire Nuclear Bulge and explains the relatively low average extinction towards the Galactic Centre. 
In addition, we find 3xl0 7 Mq of cold and dense material outside the Nuclear Bulge at positive longitudes and 
1 x 10 7 M© at negative longitudes. This material is not heated by the stars in the Nuclear Bulge and gives rise to 
the observed asymmetry in the distribution of interstellar matter in the Central Molecular Zone. 

Key words, dust, extinction - ISM: structure - Galaxy: centre - Galaxy: structure - Infrared: ISM: continuum 



1. Introduction 

Our present knowledge of the Galactic Centre (GC) 
Region has recently been reviewed by, e.g., Blitz et al. 
(1993), Genzel et al. (1994), Mezger et al. (1996; hereafter 
MDZ96), and Morris & Serabyn (1996). The gas dynamics 
of the inner Galaxy has recently been re-investigated by 
Englmaier & Gerhard (1999). To keep this introduction as 
concise as possible, we quote here only the most relevant or 
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recent papers and refer otherwise to the corresponding sec- 
tions of MDZ96 and references therein. For the same rea- 
son we use a number of abbreviations whose meanings are 
explained in Table |l|. For consistency, we adopt a distance 
to the Galactic Centre of i?o = 8.5 kpc throughout this 
paper, although recent studies suggest a somewhat lower 
value (e.g., McNamara et al. 2000: Ro = 7.9±0.3kpc; see 
also Rcid 1993). 

MDZ96 classify the centre of our Galaxy as a mildly 
active Seyfert nucleus. Although the presence of a black 
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Table 1. Abbreviations used in this paper 



Abbreviation 


Meaning 


BH 


Black Hole 


BW 


Baade's Window 


CMZ 


Central Molecular Zone 


FIR 


lar-mtrared (~30— 300 fxm) 


b WHM 


lull Width Halt Maximum 


GB 


Galactic Bulge 


GC 


Galactic Centre 


GD 


Galactic Disk 


GMC 


Giant Molecular Cloud 


HPBW 


TT IT T"Ti T~i TTT"' li 1 

Halt Power Beam Width 


ri WrlM 


Half Wiatn Hall Maximum 


ISM 


Interstellar Matter 


KLF 


K-band Luminosity Function 


MIR 


mid-infrared (~7— 30) fim) 


MS 


Main Sequence 


NB 


Nuclear Bulge 


NMD 


Nuclear Molecular Disk 


NSC 


Nuclear Stellar Cluster 


NSD 


Nuclear Stellar Disk 


NIR 


near-infrared (~1— 7/im) 


PAH 


Polycyclic Aromatic Hydrocarbon 


SED 


Spectral Energy Distribution 


ZL 


Zodiacal Light 



hole of ~ 2.6 x 10 6 Mq is strongly supported by recent 
observations (e.g., Eckart & Genzel 1998) it has also be- 
come clear that most of the activity in the centre of our 
Galaxy is due to massive star formation in the central 
parsec (MDZ96, Sect. 5). 

The mass of the central region of our Galaxy (R < 
3 kpc) is dominated by the Galactic Bulge which consists 
mainly of old, evolved stars. Both the stellar near-infrared 
(NIR) surface brightness distribution and the kinematics 
of the gas suggest that the Galactic Bulge has a bar struc- 
ture with its near end in the first galactic quadrant (i.e., at 
positive I) (e.g., Binney et al. 1991; Blitz & Spergel 1991; 
Weiland et al. 1994). As a consequence of the dynamics 
caused by the gravitational potential of this bar, the re- 
gion around the co-rotation radius Rqr ~ (3.5 ± 0.5) kpc 
is de-populated of gas, since Interstellar Matter (ISM) 
is transported efficiently inward from the Galactic Disk 
Molecular Ring at the bar's outer Lindblad resonance at 
i?OLR ~ (4 ± 0.5) kpc. Inside R < 2 kpc the gas settles 
on closed elongated (X\) orbits. This gas is observed as 
a tilted disk of atomic hydrogen (R < 1.5 kpc), usually 
referred to as the "HI Central Disk" (or "Nuclear Disk") 
(M H ~ 4x 1O 7 M ; Burton & Liszt, 1978). Inside the inner 
Lindblad resonance the X\ orbits become self-intersecting, 
and shocks and angular momentum loss compress the gas 
into molecular form and drive it further inward where 
it finally settles on more circular, stable X 2 orbits (e.g., 
Englmaier & Gerhard 1999). The existence of a distinct, 
unusually dense molecular cloud complex in the central 



few hundred pc of our Galaxy, often referred to as the 
"Central Molecular Zone", (CMZ) is well-established since 
the early 1970s (e.g., reviews by Genzel & Townes 1987 
and MDZ96). The gas distribution in the CMZ is highly 
asymmetric with most of the mass being located at posi- 
tive longitudes and positive velocities (see MDZ96, Sects. 
2.2.3 and 3.4). One of the most remarkable features in the 
I - v plane of the CMZ is the "180-pc Molecular Ring" 
which is hypothesized to be a shock region between the 
innermost stable X\ orbit of the bar and the more circu- 
lar X 2 orbits in the centre (Binney et al. 1991), but was 
also interpreted as an expanding molecular ring or shell 
(e.g., Scoville 1972; Bally et al. 1987; Sofue 1995b). 

The large and dense stellar complex in this region was 
originally thought of as the innermost part of the more 
extended Galactic Bulge, with its population of old and 
evolved stars. Various observations, e.g. ongoing star for- 
mation, the presence of ionizing stars, and its extraordi- 
nary high surface brightness suggest that this region is 
distinct from the old Galactic Bulge (as originally pro- 
posed by Serabyn & Morris 1996) and may be associated 
with the CMZ. MDZ96 therefore call the innermost region 
R < 300 pc "Nuclear Bulge" (NB), and clearly distinguish 
it from the "Galactic Bulge" (GB). Due to its relative 
proximity, the physical characteristics of the NB of our 
Galaxy can be studied in detail. However, the layer of in- 
terstellar dust in the Galactic plane restricts observations 
of the NB to wavelengths A > 2.2 /jm and the edge-on pro- 
jection makes a derivation of the true three-dimensional 
morphology difficult. 

We have begun a systematic investigation of the phys- 
ical characteristics of the NB the results of which are be- 
ing published in a series of papers. Papers I (Philipp et 
al. 1999a) and II (Mezger et al. 1999) analysed the stel- 
lar population of the central ~30pc. Based on a high- 
resolution A2.2^m survey, we determined the K-band lu- 
minosity function (KLF) of the central 30 pc and inter- 
preted it in terms of a present-day bolometric luminos- 
ity and mass function. Main Sequence (MS) stars with 
masses < 1 M Q account for ~90% of the dynamical mass, 
but only for 6% of the K-band flux density. MS stars with 
masses > 1 Mq account for ^6% of the dynamical mass 
and a similar percentage of the integrated K-band flux 
density, but are responsible for ~80% of the bolometric 
stellar luminosity as well as the ionization of the observed 
HII regions. The bulk of the K-band emission comes from 
stars evolved from the MS such as giants, supergiants, and 
Wolf-Rayet stars. We find a deficiency of low-mass stars 
within the central 1.25 pc and indications of a high star 
formation activity during the past 10 7 — 10 8 years. 

Here, in paper III, we analyse the large-scale distribu- 
tion of stars and ISM in the NB using IRAS and COBE 
DIRBE data. The size scales addressed in this paper span 
tens to hundreds of parsecs. The complex structure in the 
central few parsecs, including, e.g., the circumnuclear disk, 
the mini-spiral, and the central cavity, is not addressed. 
In a succeeding paper, based on ground-based single-dish 
mm observations, we investigate the morphology and kine- 
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matics of ISM in the central part of the NB in more detail 
(Zylka et al., in prep.). 

The paper is organized as follows: Section ^ describes 
briefly the observational data used in this paper. Section 
U describes the data reduction and analysis. In Sect, i we 
present the basic observational results, and in Sect. 5 we 
derive physical characteristics of the NB and the CMZ. 
Section ^| summarizes the results in terms of a coherent 
picture of the NB. 

2. The data 

This paper is mainly based on IRAS (Infrared 
Astronomical Satellite) and COBE (Cosmic Background 
Explorer) data covering the wavelength range from 2.2 fim 
to 240 /im, which were obtained through Sky View on the 
WorldWideWeb. Additional data at mm and cm wave- 
lengths are used for comparison. Characteristics of the 
surveys of the Galactic Centre Region used in the fol- 
lowing analysis are given in Table ^, followed by a more 
detailed description of the data. Throughout this paper 
we use Galactic coordinates (Z 11 , 6 n ) and all maps are cen- 
tered on l u ,b u = 0°,0° (and not on SgrA*, which is at 
l n ,b 11 = -0?054,-0?046). 

2.1. IRAS data 

The IRAS data include all data distributed as part of the 
IRAS Sky Survey Atlas (ISSA) and were processed by the 
Infrared Processing and Analysis Center (IPAC) to a uni- 
form standard with the Zodiacal Light (ZL) already sub- 
tracted. For the present investigation, ISSA maps in all 
four bands were used (i.e., A = 12, 25, 60, and 100 /im). 
These maps have a mean angular resolution of 2', al- 
though IPAC has optimized the processing of these data 
for features of size 5' or more. The 100 /im map is heavily 
saturated in the central part of the NB, and the 60 /im 
map may be partially saturated towards the SgrA com- 
plex. At each of the four wavelengths we obtained maps 
of 15° x 10° with a pixel scale of 1'.5, and 40° x 40° with 
a pixel scale of 5', respectively. 

2.2. COBE data 

The COBE Diffuse Infrared Background Experiment 
(DIRBE) was planned primarily to investigate the cosmic 
infrared background radiation in 10 filter bands centered 
at 1.25, 2.2, 3.5, 4.9, 25, 60, 100, 140, and 240 /im (for 
details see the COBE DIRBE Explanatory Supplement 
1997; Hauser et al. 1998; Kelsall et al. 1998). The COBE 
DIRBE maps cover the entire sky and provide an estimate 
of the infrared intensity at each pixel and wavelength band 
based on an interpolation of the observations made at var- 
ious times at solar elongations close to 90°. To first order, 
these maps depict the sky as if it were observed through a 
temporally constant interplanetary dust foreground, thus 
enabling straightforward modeling and subtraction of the 
ZL. The DIRBE instrument has a 42' x 42' instantaneous 



field of view. For this study, maps in all 10 bands were 
obtained. The maps have an angular resolution of ~0.7° 
and a pixel scale of ~0.35° (21'). 

2.3. Complementary millimetre and radio continuum 
data 

Millimetre continuum maps at 11" resolution of the cen- 
tral region of the NB were obtained with the IRAM 30-m 
telescope during several observing runs between 1995 and 
1999 (see Table ||). These maps are only used for compar- 
ison with the COBE and IRAS data and will be presented 
and discussed in detail in Paper IV by Zylka et al. (in 
preparation) . 

In addition, we obtained radio continuum maps of the 
GC region at 3, 6, 11, and 21 cm (Fig. |) from the MPIfR 
Survey archive. These data are already published else- 
where (see Table ||) and the 11cm and 21cm maps are 
already decomposed into source and background compo- 
nents. From the total intensity 3 cm and 6 cm maps, we 
subtracted the background emission to obtain the source 
contribution. Due to the different beam sizes, observ- 
ing modes, and data reduction procedures, the four data 
sets may not recover extended emission in the same way. 
Therefore, the derived flux densities given in Table | and 
shown in Fig. [l0]may contain systematic uncertainties and 
are not suitable to derive physical characteristics of the 
gas from flux density ratios. We use these maps only to 
compare the morphology of the ionized gas with that of 
the stars and dust in the NB. 



3. Data analysis 

3.1. Sources of radiation and extinction 

In order to derive surface brightness maps of the GC re- 
gion from the data, the total surface brightness maps were 
decomposed into the individual contributing sources of ra- 
diation and extinction which were, in turn, modeled and 
interpolated over the GC region before they were sub- 
tracted from the total surface brightness maps. Figure [l] 
illustrates the different Galactic features which contribute 
to emission and absorption in the direction of the GC 
and Table | summarizes the different sources dominating 
emission and absorption in different wavelength regions. 
Figure || illustrates that, particularly at NIR wavelengths, 
the resulting maps of the NB depend crucially on proper 
decomposition, extinction corrections, and modeling of the 
different galactic features. Therefore, this Section gives a 
detailed description of the data analysis process. 

The major emission sources at NIR wavelengths are 
stars in the Galactic Disk (GD) and GB. In addition to 
these features, the ZL from interplanetary dust in the so- 
lar system contributes to the extended near- and mid- 
infrared (MIR) emission. Extinction by dust in the GD 
and in the GC region itself prevents observations of the 
stellar component of the NB at wavelengths shorter than 
2 /im, and must be taken into account for wavelengths as 
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Table 2. Surveys of the GC referred to in this paper 



Wavelength range 


Angular 


Observations 


Telescopes 


References 






resolution 




(see References) 




NIR - FIR 


1.25 /im - 240 /xm 


0.7° 


10 continuum bands 


COBE DIRBE 


Hauser et al. (1991) 


MIR - FIR 


12 fim - 100 /im 


1'.5 


4 continuum bands 


IRAS 


IRAS Sky Survey Atlas (ISSA) 


Millimetre 


800 /xm 


30' 


continuum map 


CSO 10.4 m 


Lis & Carlstrom (1994) 




1.2 mm 


11" 


continuum map 


IRAM 30 m MRT 


Zylka et al. (in prep.) 




2.6 mm 


8'.8 


12 CO (1-0) map 


CSMT 1.2 m 


Bitran et al. (1997) 


Radio 


3 cm (10 GHz) 


3' 


continuum map 


Nobeyama 45 m 


Handa et al. (1987) 




6 cm (5000 MHz) 


4'.1 


continuum map 


Parkes 64 m 


Haynes et al. (1978) 




11cm (2695 MHz) 


4' 3 


continuum map 


Effelsberg 100 m 


Reich et al. (1990b) 




21cm (1408 MHz) 


<)'.4 


continuum map 


Effelsberg 100 m 


Reich et al. (1990a) 



+ 




i 1 1 1 1 1 1 1 1 ^ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

B 10 14 



Fig. 1. Schematic presentation of the different galactic features which contribute to emission and absorption in the 
direction of the Galactic Centre Region. The grey-scale image shows the average hydrogen density in the Galactic 
Disk (GD) (0.1 to 1 cm" 3 ; slice through the disk model described in Appendix |b|) . Narrow dark features correspond 
to spiral arms which have a 3 to 30 times enhanced average density. Subscripts 'fg' and 'bg' refer to the foreground 
and background parts of the GD with respect to the GC. The dashed line denotes the Galactic mid-plane and solid 
lines mark the density FWHM thickness of the gaseous disk (H 2 + HI + HII) . The interplanetary dust ring of the solar 
system (Zodiacal Light) is not shown. Note that the z scale is stretched with respect to the Rq scale. The Galactic 
Bulge (GB) is not drawn in its full extent in latitude. The Nuclear Bulge (NB) is shown on scale. 

Table 3. Sources which dominate in different wavelength regions the emission and absorption in and towards the NB 



Wavelength region 


A range 


Emission due to 




Absorption due to 


NIR 


1 jiim - 7 /xm 


stars 




dust 


MIR 


7 /xm - 30 /xm 


hot dust (VSGs) 


, PAHs 


dust 


Submm/FIR 


30 /xm - 3 mm 


"classical" dust 


grains 


dust, negligible for A > 100 /xm 


Radio 


3 mm - ~6 cm 


free- free 








> ~6 cm 


synchrotron 




free- free 



long as A ~ 100 /xm (see Table ||). Most of the extinction 
in the GD arises from the Disk Molecular Ring at galac- 
tocentric radius Rq ^4 — 5 kpc. The far sides of the GB 
and GD also suffer from extinction by dust in the GC re- 
gion. In addition to these large-scale absorption features, 
some compact Giant Molecular Clouds (GMCs) along the 
line of sight appear as discrete patches of high extinction. 



These clouds could not be modeled individually, but had 
to be masked out before the models were fitted. Since GD 
and GB have nearly exponential emission profiles, proper 
modeling was essential to obtain reliable surface bright- 
ness interpolations over the GC region. 

At wavelengths A > 7/xm, dust emission (i.e., mainly 
re-radiation of absorbed stellar emission) from the GD 
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Fig. 2. Latitude profiles of the observed total surface 
brightness at b = 0° (left panels, thick grey lines) at 
a) A2.2/im, b) A25/im, and c) A240/im as observed by 
COBE DIRBE. Black solid, dashed, and dotted lines show 
the contributions of NB, GB, ZL, and GD, respectively, to 
the observed surface brightness. Right panels d) through 
f) show the final dereddened latitude profiles of the NB 
(solid lines). In addition, the dereddened 2.2/xm surface 
brightness profile of the Galactic Bulge is shown as dashed 
line in panel d). 

and GC region dominates over stellar emission. The 
MIR emission is generally dominated by UV-excited, 
non-transiently-heated very small grains and Polycyclic 
Aromatic Hydrocarbons (PAH's) with a possible contribu- 
tion by hot dust in circumstellar shells and photospheric 
emission from cold luminous supergiants. The bulk of the 
dust (> 99% by mass) in the GD and NB is relatively 
cold (~ 15-30 K). Therefore, dust emission as tracer of the 
mass of ISM is best observed at FIR and submm wave- 
lengths. In the range 3 mm < A < 6 cm, free-free emission 
from gas ionized by early-type MS stars dominates the 
radio emission of the NB. Synchrotron radiation from rel- 
ativistic electrons becomes the strongest radio component 
for A > 6 cm, but free-free absorption has to be taken into 
account already at cm wavelengths, depending on the free- 
free emission measure. 

The total surface brightness 7 to t observed in the direc- 
tion of the GC can be written as: 



Itnt — I' 



ZL 



I, 



GD 



GB 



^nb + e~ 



('gI + 'gd)) (1) 



ixx are total "observed" intensities, i.e., integrated along 
the line-of-sight through the corresponding components, 



but not corrected for extinction by dust inside these com- 
ponents. For the meaning of the subscripts we refer to 
Table [l| In addition we use the superscripts "fg" = fore- 
ground and "bg" = background (with respect to the GC 
as seen from the Sun). Surface brightness / and optical 
depth r are functions of wavelength A, galactic longitude 
I, and latitude b. 



3.2. Goal and strategy 
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Fig. 3. Flow chart of the COBE data reduction steps 



The main goal of the data reduction was to derive sur- 
face brightness maps of the NB at different wavelengths as 
would be seen with no other emitting or absorbing mat- 
ter present in the Galaxy. The other data mentioned in 
Sect. H were reduced by its respective authors and used as 
published (references in Table ||) . 

For reasons mentioned in Sect |3.l| , surface brightness 
and extinction maps of the different galactic features had 
to be derived and modeled partially in an iterative way 
before they could be subtracted from the total surface 
brightness maps. The correct order of applying the differ- 
ent corrections to the data is described by Eq. ([!]) . Figure 
[| displays schematically the data flow with its recursive 
iterations. The general data reduction steps described in 
the following subsections are: 
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Photometric corrections for the IRAS and COBE hlter 
transmission characteristics (Sect. |3.3| ). 
Model and subtract the ZL (Sect. |3.4|). 



Model the dust distribution in the GD (Sect. 3.6 and 
Appendix [b]). 

Model the NIR surface brightness of the GB (Sect. 



5. Derive extinction maps due to dust in the GD and in 



the GC region (Sect. 3.5). 

Fit, correct for extinction, and subtract the contribu- 
tions by stars and dust in the GD (Sect. [D]). 
Correct for extinction and subtract contributions by 
stars in the GB (Sect. |3~7j) . 

Correct the final NIR maps of th e NB for extinction 
by dust inside the NB (Sect. |3.5.2| ). 
Decompose ISM in the CMZ into ISM located inside 
and outside the NB (Sect. 5.5). 



3.3. Photometric corrections 

Colour correction factors for COBE (Appendix A of the 
COBE/DIRBE Explanatory Supplement) and IRAS maps 
(Wheelock et al. 1994) were derived for the peak flux den- 
sity ratios toward the NB, but were applied to the entire 
images. The maximum correction had to be applied to 
the COBE 25 ^m band (23%). For all other bands, cor- 
rection factors are smaller then 5% of the peak flux den- 
sities. In addition, the COBE DIRBE to IRAS photomet- 
ric transformation was applied to the IRAS ISSA maps 
(COBE/DIRBE Explanatory Supplement; only gain fac- 
tors were used) in order to account for systematic photo- 
metric uncertainties in the IRAS calibration (particularly 
in the 60 and 100 /xm bands). 



3.4. Zodiacal Light subtraction 

ZL emission was already subtracted from the IRAS data 
issued by IPAC. In the COBE maps we modeled the ZL 
outside the Galactic Disk (between \b\ = 11 . . . 14°), in- 
terpolated linearly over the region of the disk, and sub- 
tracted the resulting maps from the COBE maps in the 
wavelength bands between 2.2 /im and 100 /im. 



3.5. Extinction corrections 

Two major extinction-causing features have to be consid- 
ered in the process of decomposing the observed surface 
brightness maps into individual galactic components (see 
Fig. and Eq. (§): 

1. Dust located at the near side of the GD absorbs light 
from the GB and NB, as well as from the far side of 
the GD. 

2. Dust in the CMZ absorbs light from the far sides of 
GB and GD, respectively. 

The dust opacity spectrum used thr ough out this paper is 
described in Appendix |a| (see Table A.l ). 



3.5.1. Extinction by dust in the Galactic Disk 

The foreground extinction due to dust in the GD can be 
derived in two ways: 

1 . From the observed reddening of the stellar NIR emis- 
sion from the GB, and 

2. from the optically thin FIR dust continuum emission 
of the GD. 

Neither of the two methods is straightforward. The first 
method requires the knowledge of the effective colour tem- 
perature of the GB and is very sensitive to how the GD 
emission is subtracted and where the extinction map is as- 
sumed not to be strongly affected by bright star-forming 
regions or very opaque individual GMCs. The second 
method requires a three-dimensional model of the dust 
distribution in the GD because the observed FIR contin- 
uum emission arises from lines of sight through the entire 
disk while extinction towards the GC is caused by dust on 
the near side only. 

Method 1: Under the assumption that the intrinsic NIR 
colour temperature (i.e., the intrinsic surfa ce b rightness 
ratio I1/I2) of the GB is uniform (see Sect. fT^) and can 
be derived at high latitudes where extinction by dust in 
the GD is negligible, a foreground extinction map can be 
derived from the observed NIR surface brightness ratio 



1.086- 



Ky 



K 2 - Kl 



mi 5 



J-2 



(2) 



Here, Ki are the dust opacities at the corresponding wave- 
lengths (see Table A.l). Due to the problems mentioned 
above, only an average latitude extinction profile could 
be derived rather than a complete extinction map. The 
intrinsic NIR surface brightness ratios I1/I2 of the GB 
for the A 2.2, 3.5, and 4.8 fim DIRBE bands were derived 
at \b\ > 3.5° from the disk -sub tracted (but not extinction- 
corrected) maps (see Sect. 3.6). They correspond to an av- 
erage effec tive temperature of stars in the GB of ^4000 K 
(see Sect. |4.2|) . This method yields a visual extinction 
between the Sun and the front side of the GB/NB of 
Ay ~ 14 mag. The corresponding average latitude extinc- 
tion profile is shown in Fig. ^. 

Method 2: In order to derive the extinction by dust on 
the near side of the GD from the optically thin FIR emis- 
sion arising from lines of sight through the entire disk, we 
developed a three-dimensional model of the dust distribu- 
tion in the GD, which is described in Appendix ||. Our 
best model yields Ay ~ 15 mag towards the front side of 
the NB and the latitude Ay profile at I = 0° shown in 
Fig. |. 

These two methods are completely independent. They 
yield similar, but systematically slightly different latitude 
extinction profiles. Since both methods have their weak 
points, we finally used an extinction map which was ob- 
tained by scaling the extinction map derived by method 
2 to the average of the two extinction profiles at I = 0° 
which is also shown in Fig. 0. 



Launhardt et al.: Large-scale characteristics of the Nuclear Bulge 



7 



3.5.2. Extinction by dust in the Galactic Center Region 

Since radiation from the far sides of GB and GD suffers 
from extinction by dust in the GC region, the modeled 
or interpolated surface brightness maps of these features 
have to be corrected for this extinction before they can be 
subtracted from the raw maps. Assuming a homogeneous 
and axisymmetric distribution of dust in the NB, the re- 
sults from Sect. 11 (see Fig. ||) suggest a visual extinction 
of ~ 30 mag through the entire NB. In contrast, we derive 
Ay(NB) ~ 200 mag from the FIR dust continuum emis- 
sion, assuming a homogeneous dust distribution (Sect. 5.6) 
. Uncertainties in our dust model may contribute to, but 
cannot fully explain this discrepancy. Rather we argue in 
Sect. 5.6, that the ISM in the CMZ is extremely clumpy 



with most of the mass being concentrated in small, ultra- 
opaque molecular clouds which cover ~ 10% of the area of 
the NB and block the stellar radiation nearly completely. 
In Sect. 4.3.2 we derive Ay(NB) ~ 20 mag for the average 



diffuse extinction through the NB. Therefore, the far side 
parts of GB and GD were corrected with an extinction 
map which was derived by scaling the total hydrogen col- 
umn density map of the CMZ (Fig. p~6|b) to Ay = 20 mag 
plus a 10% blockage factor within the A^h — 2x 10 22 cm~ 2 
contour. The assumption that the degree of dumpiness of 
ISM is the same inside and outside the NB may not hold 
true completely, but is certainly a better approximation of 
the physical conditions than assuming homogeneous mat- 
ter distribution. 

To derive the size and morphology of the stellar NB, 
the NIR emission from the NB itself was, in addition, cor- 
rected for extinction by dust inside the NB. In Sect. |5.5| 
we decompose the total hydrogen column density map of 
the CMZ into material located inside and outside the NB. 
An intrinsic extinction map of the NB was then derived 
by scaling its A^h map (Fig. |l6|d) in the same way as above 
(i.e., v4y(NB) = 20 mag plus 10% blockage). In correcting 
the NIR surface brightness of the NB for intrinsic extinc- 
tion, we assumed that stars and diffuse dust within the 
NB are homogeneously mixed. Intrinsic (/*) and observed 
surface brightness (I' v ) are then related by: 



I' v {\) = 0.9 (1 



, 4* (A) 



(3) 



3. 6. Galactic Disk subtraction 



The "observed" surface brightness distribution of the GD 
at all wavelengths was derived by interpolating the to- 
tal observed surface brightness between selected regions 
in the longitude intervals 3.8°< 1 1 1 < 20° (i.e. within the 
inner tangent to the Disk Molecular Ring, but outside the 
NB). The regions used for deriving the average disk emis- 
sion profile were determined by comparing NIR and FIR 
maps in order to avoid areas which deviate strongly from 
the average disk profile, like heavily FIR-emitting or NIR- 
absorbing GMCs, bright star-forming regions, or gaps in 
the dust disk. At NIR wavelengths (DIRBE bands 1 to 
4), GD and GB are superimposed at longitudes \l\ < 20° 



and the surface brightness of the GD cannot be derived 
independently. Therefore, it was first modeled further out 
(\l\ > 20°) and subtracted from the total emission before 
modeling and subtracting iteratively the surface bright- 
ness distributions of GB (see Sect. 3.7) and inner GD. 
Before the GD emission was finally subtracted from the 
total surface brightness maps, it was decomposed into con- 
tributions from the near and far sides (considering extinc- 
tion by dust in the near side of the GD) and corrected for 



extinction in the GC region as described in Sect. 3.5 



3. 7. Galactic Bulge subtraction 

To derive NIR surface brightness maps of the NB, all 
contributions from the GB must be subtracted. Here, we 
assume that the GB and NB are morphologically differ- 
ent features and that the central surface brightness pro- 
file of the GB can be derived from its outer parts. Since 
the surface brightness profile of the GB has a centrally 
peaked exponential (or similar) shape and is not sym- 
metric, proper extinction correction and interpolation over 
the central region is crucial. When properly corrected for 
foreground extinction, the often-referred peanut shape of 
the surface brightness distribution nearly disappears and a 
box-shaped asymmetric ellipsoid remains (Fig. ^|a). Since 
there is striking evidence that the GB is a triaxial stellar 
bar with its near end at positive longitudes (e.g., Binney 
et al. 1991; Blitz & Spergel 1991; Weiland et al. 1994), 
its surface brightness distribution cannot be modeled by a 
simple symmetric function. Therefore, we followed for the 
most part the approach of Dwek et al. (1995), Stanek et 
al. (1997), and Freudenreich (1998), and modeled the GB 
as a triaxial ellipsoidal bar with a centrally peaked volume 
emissivity p. We have tried sech 2 , Gaussian, and exponen- 
tial volume emissivity profiles and found the exponential 
models, represented by: 



p oc exp (—Rs) 



(4) 



to give the best fits. The shape of the bar was chosen 
to be a "generalized ellipsoid" (Althanassoula et al. 1990; 
Freudenreich 1998) represented by the effective radius i? s : 



Re 



R 



\Z'\ 



with 



\X'\ 



(5) 



(6) 



where a x , a y , and a z are the scale lengths and C± and Cm 
the face-on and edge-on shape parameters. The orienta- 
tion of the bar was characterized by an in-plane tilt 4> x (an- 
gle between bar's major axis and the line Sun - Galactic 
Centre) and an out-of-plane tilt 9 Z . To derive the surface 
brightness, the models were integrated as seen from the 
position of the Sun which was adopted to be Z$ — 16 pc 
above the galactic mid-plane (e.g., Freudenreich 1998). 

In the GD-subtracted COBE NIR maps of the GB 
corrected for extinction by dust in the GD (Sect. 3.5), 



8 



Launhardt et al.: Large-scale characteristics of the Nuclear Bulge 



both the region of the CMZ and of individual opaque 
GMCs were masked out before the models were fitted to 
the observed surface brightness distribution. Parameters 
of our best-fit GB model are given in Table ||. The Spectral 
Energy Distribution (SED) of the GB is derived and anal- 
ysed in Sect. 4.2 (Fig. ^). Figure || compares the model 
with the observed surface brightness distribution. 

The modeled surface brightness maps of the GB were 
then corrected for extinction by dust in the GC region 
(i.e., the CMZ): 



I'r 



GB — i GB ' 



(1 + e - TCMZ ) 



(7) 



(see Sect. 3.5.2 for tcmz) and subtracted from the COBE 
maps. 

4. Basic Parameters derived from the COBE 
DIRBE and IRAS ISSA maps 

4.1. Extinction towards the Galactic Centre Region 



profile, which we finally used, can be approximated by 
the sum of two sech 2 functions: 



A v (b) 



i=l 



A Vt (0)seclr — 



(8) 



with Avi(O) ~ 11 mag, a bl ~ 0.78°, A V2 (0) ~ 3.5 mag, 
and (Tf,2 — 3.3°. The total mid-plane extinction to- 
wards the front-sides of GB and NB is Ay = 14.5 ± 
2 mag. Compared to the well-established value Ay(Sun- 
GC) ~ 30 mag (e.g., Lebofsky & Rieke 1987; Catchpole et 
al. 1990), our result suggests that ^50% of the absorbing 
dust column between the Sun and the GC is located in the 
GD and the other 50% is located in the GC region. Figure 
[I| also shows the extinction profile derived by Catchpole et 
al. (1990) from colour-magnitude diagrams of stars in the 
direction of the GC. This profile represents the extinction 
to stars in the immediate vicinity of the GC and accounts 
for extinction by dust in the GD and in the GC region. 

4.2. Galactic Bulge 



30 



20 



Catchpole et al. (1990) x ^» 



used for 
dercdccTiinc 



GB reddening (NIR) 
( m et hod 1 ) 




Fig. 4. Visual extinction towards the front-side of 
Galactic and Nuclear Bulge, due to dust in the Galactic 
Disk. Dotted curve: Extinction profile derived from 3D 



model of the dust distribution in the GD (Sect. 3.5.1 



Method 2). Dashed curve: Extinction profile derived from 



reddening of the GB (stellar) NIR emission (Sect. 3.5.1 
Method 1). Solid curve: Average of the two former pro- 
files, which we use to correct the NIR maps of the NB. 
Grey dots: Extinction profile derived by Catchpole et al. 
(1990) from colour-magnitude diagrams of stars in the di- 
rection of the GC. This profile relates to the Centre of the 
Galaxy, i.e. it includes extinction by dust inside the NB. 

In Fig. ||, the latitude profiles of extinction between 
the Sun and the fro nt sides of GB and NB at I = 0° , as 
derived in Sect. |3.5|, are shown. The average extinction 



Table 4. Galactic Bulge model parameter values 
(see Eqs. (@) to (§)) 



Parameter 


Value 


Sun dist. from Galactic Center Ro a ' 


8.5 kpc 


Sun dist. from Galactic Plane Zo 


16 pc 


Bar X scale length a x 


1.1 kpc 


Bar Y scale length a y 


0.36 kpc 


Bar Z scale length a z 


0.22 kpc 


Bar face-on shape parameter C± 


1.6 


Bar edge-on shape parameter C\\ 


3.2 


Bar in-plane tilt angle cj> m 


15° 


Bar out-of-plane tilt angle Z 


-1° 


&, (1.25 pm) 


(0.92±0.32) MJy 


S v (2.2 pm) 


(0.82±0.18) MJy 


S„(3.5 fim) 


(0.42±0.07) MJy 


(4.9 pm) 


(0.21±0.02) MJy 


Average effective temperature T e g 


(4400±400) K 


Total Luminosity Lboi 6 ' 


(1.0±0.3)xl0 10 L 



a ' Fixed input parameter 

6) Integration of the SED (Fig. 



The parameters of our 'best-fit' 
the GB (Sect. O and Eqs. (0) to (S 



triaxial bar model of 
J)) are summarized in 
Table |^, and the result is shown in Fig. ||. Our model yields 
similar axis ratios, (a x : a y : a z ) ~ (3 : 1 : 0.6), and projected 
surface brightness profiles as the best-fit models of Dwek 
et al. (1995) and Freudenreich (1998). The tilt angle of the 
bar's major axis against the line of sight in our model is 
4> x ~15°. Other studies favor a somewhat larger tilt angle 
of 16-25° (e.g., Englmaier & Gerhard 1999). Note that 
our goal was to obtain a good fit to the projected surface 
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5 0-5 
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Fig. 5. 2.2 /im surface brightness distribution of the 
Galactic and Nuclear Bulge. In order to achieve the best 
signal-to-noise ratio, a weighted average of the 2.2, 3.5, 
and 4.9 /im maps, scaled to the 2.2 /im surface brightness, 
is shown. Contributions from the Galactic Disk are sub- 
tracted and the emission is dereddened for extinction by 
foreground dust, a) Contour maps of the observed dered- 
dened (thick grey lines) and modeled (dashed lines) NIR 
surface brightness distribution. Levels are at 5, 10, 20, 
60% of modeled peak surface brightness, b) Latitude pro- 
file at I = 0° of the observed dereddened (grey dots) and 
modeled (dashed line) NIR surface brightness distribution 
of the GB. 



brightness distribution of the GB in order to subtract it 
from the COBE maps, rather than to derive a detailed 
model of the three-dimensional morphology of the bar. 
We did not explore the whole parameter space and we do 
not claim that our model is unique. 

In Figure ^, the integrated SED of the GB is shown. 
The NIR flux densities were derived by integrating the 
models at 2.2, 3.5, and 4.9 /im within I — ±20° and b = 
±10° and are listed in Table ||. A weighted least-square 
black-body fit to these points yields an average effective 
temperature of the Bulge stars of T cff (GB) = 4400±400K. 
Since we found no evidence for a colour gradient in the GB, 
a lower limit to the average effective stellar temperature 
can be derived from the uncorrected NIR surface bright- 
ness ratios at high latitudes where the extinction is low 
(e.g., Baade's Window: l,b -l°,-3.9°, A K ~ 0.13 mag). 
The best black-body fit to these points yields ~ 3850 K 
which is in good agreement with our estimate of T e g (GB). 
We tried different bar models and fitting routines; the de- 




v [Hz] 



Fig. 6. Spectral Energy Distribution of the Galactic 
Bulge. 

Black squares refer to the integrated flux density of the 
model fit (left scale) . Error bars refer to the uncertainties 
of the extinction corrections (a (Ay) ~ ±2 mag). Curve 
a) shows the best black-body fit to these points weighted 
with the uncertainties (T c g ~ 4400 K). Filled circles re- 
fer to the surface brightness at Baade's Window (BW) in 
the disk-subtracted, not extinction-corrected COBE maps 
(right scale). Curve b) shows the best black-body fit to 
these points (T oS ~ 3850 K.) 



rived total flux densities were nearly independent of the 
particular model. 

Our result agrees well with the volume emissivity ra- 
tios derived by Freudenreich (1998) for his models S and 
E. The total luminosity of the GB derived from our SED 
fit to the NIR flux densities is L GB = 1.0 ± 0.3 x 10 10 L . 
Since there is no evidence for ongoing star formation and 
the presence of hot massive stars in the GB, which would 
contribute considerable luminosity at shorter wavelengths, 
this value represents the bolometric luminosity of the GB. 
Our estimate is intermediate to the GB luminosities of 
5.3 ± 1.6 x 10 9 L o derived by Dwek et al. (1995) from 
COBE NIR observations and the 2.2 /im luminosity func- 
tion of Bulge stars and of ~ 2 x 10 10 L derived by Maihara 
et al. (1978) from their 2.4 /xm observations of the GC re- 
gion. Dwek et al. (1995) derive a total stellar mass of the 
GB ofM GB ~ 1.3±O.5xlO lo M and a mass-to-luminosity 
ratio of ~ 2 M /L . The central mass and luminosity vol- 
ume densities of the GB are pm ~ 8 ± 2 M pc -3 and 
Pl ~ 4 ± lL pc -3 , respectively. The numbers depend 
only weekly on the value of the total luminosity, which is 
more sensitive to the exact integration area and how the 
outer GB is modeled. 



4.3. Nuclear Bulge and Central Molecular Zone 

4.3.1. NIR through Radio images of the Nuclear Bulge 

Images and surface brightness profiles of the NB ranging 
from A2.2^m to 240 ^m derived from the COBE DIRBE 
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Fig. 7. Surface brightness maps of the Nuclear Bulge 
at 9 wavelengths between 2.2 /xm and 240 (xm as seen by 
COBE. ZL, emission from GD and GB, and a point-like 
NIR source at l,b ~ — 1.4°, 0° (dotted curve in the longi- 
tude profiles; see text) have been subtracted. The data are 
corrected for foreground extinction by dust in the GD, but 
not for extinction by dust inside the NB. Lowest contour 
levels are at 10% of the maximum in the NIR maps and 
5% in all other maps. Small boxes in the lower left of the 
maps show the DIRBE beam (HPBW 0.7°) at the corre- 
sponding wavelength bands. The middle and right panels 
show the corresponding longitude and latitude profiles at 
6 = 0° and 1 = 0°, respectively. 



and IRAS ISSA maps and processed as described in Sect. 
U are shown in Figs. [?] and [|. Note that these images are 
not corrected for extinction by dust inside the NB and 
that the FIR maps contain emis sion from dust in the NB 
and in the outer CMZ (see Sect. 5J3). For comparison, we 
also show radio continuum maps of the NB which were 
obtained from different data bases (Fig. ^ see Table || for 



2 a -z 
1 [deg] 

Fig. 8. Thermal dust emission from the Nuclear Bulge 
as seen by IRAS (ISSA maps, angular resolution ~2'). 
ZL and GD are subtracted and the data are corrected for 
foreground extinction by dust in the GD. Middle and right 
panels show the corresponding longitude and latitude pro- 
files as in Fig. 0. 




1 [deg] 



1D1 

b [deg] 



Fig. 9. Radio continuum emission from the Nuclear 
Bulge. These data were obtained from different data bases 
(see Table|| f° r references). FWHM beam sizes are shown 
as black circles in the lower left corners of the maps. 
Middle and right panels show the corresponding longitude 
and latitude profiles as in Fig. 0. 



references). At NIR wavelengths, the NB emits stellar con- 
tinuum radiation, mainly from red giants and supergiants. 
The NB also emits strong MIR/FIR/submm dust emission 
together with free- free and synchrotron radio emission. In 
addition, the molecular gas exhibits strong line emission 
(see Fig. [Tqa). Pertinent fit and other parameters are given 
in Table § 

The COBE NIR images show a relatively compact 
source elongated ~2-3° in longitude and unresolved in 
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latitude (upper three panels in Fig. |^). The 4.9 /Km image, 
which is least affected by extinction, shows this elongated 
structure most clearly. The 2.2 /im image exhibits some 
low-level extended emission, which may be due to uncer- 
tainties in the extinction correction and de-composition 
process and will, therefore, not be discussed here. At 
A1.25/im (DIRBE band 1), the extinction towards the 
NB and the uncertainties in the data reduction are too 
large to produce a meaningful result; therefore we don't 
show nor use the 1.25 (im image here. Visible especially 
at 2.2 /im is a source at I ~ —1.4° which has a "bluer" 
SED than that of the average NB. This source, referred 
to as IRAS 17393-3004, is the brightest NIR source seen 
in the direction of the NB. It has been investigated by 
Philipp et al. (1999b) and is classified as an M4 super- 
giant surrounded by a dust shell and located at a dis- 
tance of <4.7kpc from the sun. Therefore, we modeled 
this source by a Gaussian and subtracted it from the fi- 
nal NIR images (dotted curve in the longitude profiles in 
Fig. 0). The surface brightness longitude profiles of the 
resulting images are asymmetric and peak at I ~ +0.2°. 
Although this is marginal compared to the 0.7° DIRBE 
beam, this asymmetry is seen at all wavelengths and only 
at A > 140 /im does the peak shift towards larger /. 

The COBE MIR images are less extended in longitude 
than the NIR and FIR images, indicating that hot dust 
emission arises mainly from the inner part of the NB. The 
12 /im image shows an extended, low-level halo. Since the 
strong PAH feature at 11.3 /im lies in this band, this may 
indicate that UV-excited PAHs in and around the NB 
are more extended than "normal" warm and hot dust. 
The 12 /jm halo is also seen in the IRAS image which has 
a much higher angular resolution (Fig. ||). The IRAS 12 
and 25 /im images, too, indicate that the bulk of the MIR 
emission arises from a very compact region which includes 
the SgrA Radio Complex and the Radio Arc and Bridge 
(MDZ96; Sofue 1994; Reich 1994). 

The COBE FIR images show the same general mor- 
phology as the NIR images, i.e., the emission is extended 
in I with an asymmetry towards positive I, and is basi- 
cally unresolved in b. While at negative longitudes the 
emission extends out to I ~ —1.7° at all wavelengths, its 
extent towards positive longitudes increases with increas- 
ing wavelength up to nearly 4° at 240 /im. This indicates a 
large excess of cold dust at positive I which extends much 
further out then the stellar NIR emission and the emis- 
sion of warm dust. The excess of cold dust emission at 
I > +1.7° can also be seen in the 100 /im IRAS image 
(Fig. ||). The COBE 140 and 240 /im images likewise show 
some extended low-level emission from cold dust outside 
the NB at —2° > I > —5°. At much higher angular res- 
olution than the COBE maps, the IRAS 60 and 100 /im 
images show a very distinct and narrow ridge of emission 
between I ~ ±0.7° (Fig. ||), which is unresolved (in lat- 
itude) in the COBE maps (Fig. Q). The morphology of 
the extended low-level FIR emission in the IRAS maps 
compares well to that of the COBE images. 



Although the morphology of the radio continuum maps 
is mainly determined by emission from distinct compact 
sources and source complexes, their overall shape is sim- 
ilar to that of the MIR IRAS maps, thus indicating that 
free-free emission from ionized gas (as observed at 6 cm) 
and IR emission from hot dust have a similar distribution. 
The 11 and 21cm images are dominated by non-thermal 
synchrotron emission. 

4.3.2. The Spectral Energy Distribution of the Nuclear 
Bulge 
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Fig. 10. Spectral energy distribution of a) the cen- 
tral 120 pc and b) the entire NB. The different symbols 
represent COBE (squares), IRAS (circles), ground-based 
submm ( triangles), and radio (pentagons) data. The 
IRAS and ground-based radio data were smoothed to the 
COBE resolution. The solid lines show the best model fit 
to the observed dust and stellar continuum data (corrected 
for foreground extinction only). The dotted curves show 
the underlying dust and stellar continuum. The contribu- 
tion of radio continuum (free-free and synchrotron radia- 
tion) emission is also marked by do tted lines. Fit parame- 
ters are given in Table ||. Sect. T3^ describes how the SED 
fits were derived. Note that these SEDs do not account for 
the contribution of hot massive stars at A < 2 /im. 



To compile the radio through NIR SEDs of the central 
120 pc and of the entire NB, the corresponding central 
surface brightnesses and integrated flux densities were de- 
rived at all wavelengths from the final maps presented in 
Figs. through |[ In these maps, all contributions but the 
NB were removed and the remaining emission was cor- 
rected for extinction by foreground dust, but not by dust 
inside the NB. The total flux densities 5* ot (A) were m- 
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Table 5. Flux densities and source sizes derived from IR and radio images of the Nuclear Bulge 
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Beam-averaged surface brightness 

First valu e: cor rected for foreground extinction, (second value): additionally corrected for extinction by dust inside the NB 
(see Sect. |3.5.2| ) 

Smoothed to the COBE DIRBE resolution (0.7°) 

Emission is dominated by strong point sources. Size derived after subtraction of point sources 
The IRAS 100 /-im emission is saturated in the central region of the NB 
The submm continuum maps cover only the inner ~ 1.4° x 0.4°. 



tegrated within a polygon which completely includes the 
3% contour of the column density map (Fig. |l6lb) and 
which approximately covers the range — 2.5° < I <3.8°, 
and — 1.3° < b < 1.1°. This area partially includes dust 
in the CMZ located outside the actual NB (see Sect. |5.5| ). 
In addition, we derived surface brightnesses from the two 
submm continuum maps listed in Table |^. The IRAS, 
submm, and radio continuum maps were convolved with 
the DIRBE beam to obtain comparable surface bright- 
nesses. The central 120 pc refer to the square-shaped 0.7° 
DIRBE beam (see Fig. |7|) with an equivalent circular aper- 
ture of 0.395°, which corresponds to i? equ iv. ~ 60 pc. The 
resulting surface brightnesses and flux densities are com- 
piled in Table || and were used to construct the corre- 
sponding SEDs shown in Fig. [l(]. The error bars cor- 
respond to the usually adopted calibration uncertainties 
of -20% for the COBE and IRAS 12 to 240 ^m data 
and ~ 30% for the submm data. The uncertainties of the 
COBE NIR data were estimated from the variations when 
using different GB and extinction models (40% at 2.2 /im, 
30% at 3.5 and 4.9 /im). 

To fit the observed SEDs and derive physical param- 
eters of stars and dust in the NB, we use a simple (not 



self-consistent) radiative transfer model consisting of three 
single-temperature dust components and one stellar com- 
ponent. All four components are homogeneously mixed. 
Hot massive MS stars are not considered since their con- 
tribution to the K-band flux density is marginal (except 
for the central few pc; see Sects. 5.1 through |5~3|) and we 



don't have data points shortward of 2.2 /im. The cold and 
warm dust components represent large "classical" grains 
with their temperature distribution. Hot dust emission 
represents very small stochastically heated grains with a 
possible contribution by PAHs as well as dust in hot cir- 
cumstellar shells. The NIR part of the SED was calculated 
using the model for internal extinction in the NB described 
in Sect. 3.5.2 . The final SEDs are \ 2 fits to the observed 
flux densities with the three dust temperatures and masses 
and the effective stellar NIR temperature Tnir being free 
parameters. Here, Tnir is the effective black-body tem- 
perature derived from a fit to the completely extinction- 
corrected NIR flux densities (2.2 to 4.9 fim). Although the 
2.2 //m flux density may be slightly affected by emission 
from hot stars for which the SEDs shown in Fig. [l(] do not 
account, Tnir should approximately represent the average 
effective stellar temperature of cool evolved stars, which 
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Region (diameter): Central 1.25 pc a) Central 120 pc 6) Entire NB c) 





T 




L 


T 


M H 


L 


T 


M H 


L 




[K] 


[M ] 


[Le] 


[K] 


[M ] 


[Lq] 


[K] 


[M ] 


[Le] 


Cold dust 


40 


3.3E+2 


1.7E+5 


24.7 


3.6E+6 


5.6E+7 


21.2 


4.7E+7 


2.7E+8 


Warm dust 


150 


5.8E+0 


3.3E+6 


59 


1.8E+4 


3.6E+7 


49 


4.2E+5 


2.1E+8 


Hot dust 


350 


5.0E-3 


2.1E+5 


200 


1.1E+1 


1.1E+7 


220 


1.6E+2 


7.4E+7 


Dust (total) 


42 


3.3E+2 


3.7E+6 


24.9 


3.6E+6 


1.0E+8 


21.4 


4.7E+7 


5.5E+8 


Cool stars(obs) d)e) 












3.2E+7 






1.2E+8 


Dust + cool stars (obs) d)e) 






6.8E+6 






1.3E+8 






6.7E+8 


Cool stars(intrinsic) d ^ 


4 000 




3.1E+6 9 ' 


4600 




1.9E+8 


4400 




6.8E+8 



See Paper I, Table lb 

Corresponds to one 0?7 COBE DIRBE pixel 
' Includes dust in the immediate environment of the NB, but not the entire CMZ (— 2?5 < I < 3?8) 
Stellar NIR luminosity Z/nir (see definition below), does not account for hot (i.e., young massive) stars 
Corrected for foreground extinction, but not for extinction by dust in side t he NB 
Additionally corrected for extinction by dust inside the NB (see Sect, 3.5.2) 
Values from Paper I, Table 1. Note that Ay = 31 mag was used as total foreground extinction 
Not derived since no data points shortward of 2.2 /im available 



dominate the NIR luminosity of the NB (see Sect. 5.1). 
The contribution of hot stars in different regions of the 
NB will be further discussed in Sects. 5.1 through 5.3 



While the fit parameters for the dust emission sec- 
tions of the SEDs are nearly independent of the partic- 
ular extinction model of the NB, the results for the stellar 
(NIR) section depend strongly on the assumed distribu- 
tion of extinction within the NB. Since the populations of 
evolved stars in the Nucle ar a nd Galactic Bulge are found 
to be similar (see. Sect. 5T), we assumed Tnir(NB) ~ 
Tnir(GB). This assumption could be fulfilled by se tting 
A V (NB) ~ 20 mag (plus 10% blockage; see Sects. |3.5.2 



and |5.6| ). Higher values for the diffuse average extinction 
through the NB lead to unreasonably high temperatures 
and intrinsic luminosities (e.g., Ay(NB) = 30 mag yields 
?nir ~ 8000K and L*,nir(NB) - 3.5 x 1O 9 L ). The im- 
plication of this low diffuse extinction will be discussed in 
Sect. |5.6| . Our model yields a total intrinsic stellar NIR 
luminosity^ of i* 1 NiR(NB) ~ 7 x 10 s L , which does not 
account for the contribution from hot stars. 

The corresponding SED fits to the central 120 pc and 
the entire NB are shown in Fig. 



10 



and the results of the 
spectral decomposition and the derived masses and lumi- 
nosities are given in Table ^. For the NB we used a relative 
metallicity Z/Z & — 2 (Appendix |b], for the dust model 
used see Appendix [A|) . Note that the dust emission in this 
model, and hence the derived hydrogen masses, account 
for all interstellar matter in the central ~6°. The uncer- 
tainties and the implications of the parameter fit to the 



1 £nir is the luminosity derived from a black-body fit to the 
extinction-corrected 2.2 to 4.9 /im flux densities, assuming that 
the 2.2 /im flux is not contaminated by emission from hot stars. 



SED will be discussed in Sect. ||. The corresponding SED 
for the central 1.25 pc and fit parameters of its spectral de- 
composition given in Table || are taken from Paper I. Note 
that in Paper I Ay = 31 mag was used as total foreground 
extinction, assuming that all contributing stars are actu- 
ally located in the central parsec (no line-of-sight spread) . 
In Sect. H we make a further attempt to decompose the 
emission from the NB into contributions from different 
sub-components. 

4.3.3. Column density and mass of Interstellar Matter 
in the Central Molecular Zone 

A hydrogen column density map of the CMZ (i.e., all 
molecular material in the central kpc) can be derived from 
the optically thin 240 /im dust emission (surface brightness 
I u ) if the dust temperature Td is known: 



1 





B„(A,T d ) tohk(A) M d 



(9) 



Here, B V (X, Td) denotes the Planck function, Td the mass- 
averaged dust temperature^, mn the mass of an Hydrogen 
atom, k the mass absorption coefficient (dust opacity), 
and Mn/Md the hydrogen-to-dust mass ratio (for the dust 
model used see Appendix ^). 

Since the dust grains in the NB have a broad tem- 
perature distribution, a simple black-body colour temper- 



2 Mass-averaged dust temperature: average temperature of 
all dust components (mixed or separated along the line-of- 
sight) weighted with their respective masses. It is not neces- 
sarily identical with the effective colour temperature derived 
from flux density ratios. 
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ature derived from two FIR maps is not a good repre- 
sentation of the true dust temperature. Instead, we used 
an empirical colour temperature map derived from the 
COBE 60-to-140/mr flux ratio map of the NB. The 140 
and 60 /im fluxes are dominated by emission from cold 
(15-30K) and warm (40-100K) dust ("classical" grains), 
respectively. Hot, stochastically heated very small grains 
mainly contribute at 12 and 25 /im. The 60-to-140/im sur- 
face brightness ratio was scaled to the mass-averaged dust 
temperatures derived from the SED model fits to the cen- 
tral 120 pc and the entire NB (Fig. |l0| and Table |). We 
then assumed proportionality between Td and Sea/Sua 
in the temperature range 15. ..30K. The uncertainty of 
this method is smaller than ±1K, which results in a 30% 
uncertainty in TVh at 16 K and 10% at 25 K. Where the 
surface brightness in the 60 and 140 /im COBE maps was 
below a certain cut-off level and no reliable dust temper- 
ature could be derived, the temperature was set to 16 K, 
the lowest value derived in the map. 

The resulting dust temperature profile is shown in Fig. 
|l5| b and will be discussed in Sect. |5.5| . Mass-averaged dust 
temperatures in the CMZ are in the range 16 to 25 K. 
The corresponding dust temperature map was then used 
to convert the COBE 240 /im map (Fig. |^) into a hydrogen 
column density map (Fig. |l6|b). This map shows a convinc- 
ing similarity to the 12 CO(1-0) map of Bitran et al. (1997; 
see Fig. fl6^i), indicating that both optically thin FIR dust 
emission and 12 CO emission trace, at least on large scales, 
the same ISM. Beam-averaged column densities range 
from 1 to 15 x 10 22 cm" 2 with 7V H (0 o ,0°) - 7 x 10 22 cm" 2 
(Fig. |l5|a). The asymmetry in the dust temperature and 
column density profiles and the relation of ISM to the stel- 
lar mass distribution in the NB will be discussed in Sect. 
5J5. 

Integrating the hydrogen column density map (Fig. 
|l6|b) yields a total hydrogen mass of Mh ~ 6.8 x 10 7 M Q in 
the central 8° (~ 1 kpc) of the Galaxy. The area —2.5° < 
I <3.8° contains Mh ~ 5.3 x 10 7 M© which is in good 
agreement with the mass of 4.2 x 10 7 M Q derived from the 



SED model fit to the i nte grated flux densities (Sect. 4.3.2 
and Table ^). In Sect. 5J5 we decompose the iVn map into 



components inside and outside the NB and derive the cor- 
responding hydrogen masses and column densities. 

5. Physical characteristics of the Nuclear Bulge 

Based on the observational results presented in the previ- 
ous Section, we derive and further discuss in this section 
physical characteristics of stars and ISM in the NB and 
its immediate environment. 

5.1. The stellar population of the Nuclear Bulge 

A comparison of the model K-band luminosity function 
(KLF) of the central 30 pc with the KLF observed in 
Baade's Window indicates that the overall populations 
of low and intermediate-mass MS stars in Nuclear and 
Galactic Bulge are similar, but that the central 30 pc 



have an overabundance of K-luminous giants. These stars, 
which are interpreted as remnants of high star-formation 
activity some 10 7 to 10 8 yrs ago (e.g., Genzel et al. 1994), 
are more concentrated towards the centre than low-mass 
MS stars (Paper I, Fig. 6). Unique for the NB is a popu- 
lation of massive, young MS stars and supergiants which 
are also strongly concentrated towards the centre (Paper 
II; Genzel et al. 1996). The contribution of these stars 
to both dynamic mass and NIR luminosity is negligible 
(~ 7% and ~ 6%, respectively), but they dominate the 
bolometric luminosity in the central 30 pc (~ 80%) and 
are responsible fo r th e ionization of the thin intercloud 
medium (see Sect. 5J3). The main contributor to the stel- 
lar mass are low-mass MS stars (> 90%), which account, 
however, for only ~ 6% of the NIR luminosity. The contri- 
bution of evolved stars to the stellar mass is negligible, but 
they dominate the NIR luminosity (~ 88%; for this and 
the above quantitative estimates see Paper II). Therefore, 
an estimate of the total stellar mass from the NIR lumi- 
nosity is only possible when certain assumptions about the 
mixing of low-mass MS stars and evolved stars are made 
(see Sects. 5.3 and FO). 



5.2. Morphology of the stellar Nuclear Bulge and the 
large-scale distribution of stars 
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Fig. 11. a) Contour map of the final (observed) 4.9 /im 
map of the NB, corrected for extinction by dust in front of 
and inside the NB (thick grey contours, 10 to 85 in steps 
of 15% of the maximum surface brightness). Overlayed 
are the surface brightness contours of the best-fit model 
(shown in panel b) convolved with the DIRBE beam (thin 
black contours, same levels), b) Surface brightness map 
of the 'best fit' NIR model of the NB (see Sect. O). 



The size and morphology of the stellar NB can best 
be derived from the 4.9 /Ltm image where stellar emission 
dominates and extinction effects, i.e., uncertainties intro- 
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Table 7. Characteristics of Nuclear Stellar Cluster and Nuclear Stellar Disk 



15 













NSC: 


po (Eq. ©)W 
central lpc 3(i>) 
central lpc 2(c) 
entire NSC 


3.3 x lO 6 M /pc~ 3 
7.0 x 10 5 M Q (<i) 
1.5 x 10 6 M Q 
3 ± 1.5 x 1O 7 M 


6.6 x lO 6 L /p C - 3 
1.4 x 1O 6 L 
3.1 x lO 6 L (e) 
6±3 x 1O 7 L 


2.0 x 10 8 L e /pc~ 3 
4.2 x 1O 7 L 
9.2 x 10 7 L Q (e) 
1.8 ±0.9 x 1O 9 L 


NSD: 


inner R < 120 pc 
entire NSD 


8±2 x 10 8 M Q (/) 
1.4 ±0.6 x lO 9 M (/) 


4± 1 x 1O 8 L (9) 
6.5 ±2 x 1O 8 L (9) 


4 ± 1 x 10 s L 
6.5 ±2 x 10 s L 


total NB: 




1.4 ±0.6 x 1O 9 M 


7±2 x 10 8 L Q (h) 


2.5 ±1 X 1O 9 L 



With n = 2 for R < 6 pc, n — 3 for R > 6 pc, and an outer cut-off radius i? ou t = 200 pc 
Integrated in a sphere with R — 0.625 pc (volume V = 1 pc 3 ) 
^ Integrated in a zylinder with R = 0.625 pc throughout the entire NSC (pencil beam with area A= 1 pc 2 ) 
Genzel et al. (1997); fixed input parameter 
Paper I, Table lb; fixed input parameter 



(d) 
(<-•) 
(/) 

(s) 



(/,.) 



Derived from Lnir with M*/Lnir = 2M /L (see. Sect. 5.4) 

Derived from a two-component volume emissivity model fit of the form of Eq. (|ll|) to the extinction-corrected 4.9 /mi COBE 
map of the NB (Fig. |l|). 
Table ^| line 7; fixed input parameter 



^ Ltot includes the contribution from hot massive MS stars (see Sect. 5.3). 



duced in the data reduction and extinction correction pro- 
cess, are relatively small. Figure 11 a, shows the final 4.9 /im 
COBE image of the NB corrected for foreground extinc- 
tion and for extinction by dust inside the NB (see Sect. 
3.5.2). Corresponding longitude and latitude profiles are 
shown in Figs. |l2] e and f. Despite some sub-beam-scale 
asymmetry at the centre, the overall NIR emission from 
the NB is symmetric with respect to the position of Sgr A*. 
Since the coarse angular resolution of the COBE data pro- 
hibits interpretation of asymmetries on such a small scale, 
we assume that the large-scale distribution of stars in the 
NB is symmetric with respect to the GC. 

To derive the physical size of the stellar NB, models of 
the volume emissivity distribution were constructed, the 
resulting surface brightness distribution calculated, con- 
volved with the COBE beam, and then compared to the 
extinction-corrected 4.9 /im map. 

It was found that the central few parsec, which are 
not resolved by COBE, are dominated by the luminous 
Nuclear Stellar Cluster with an R~ 2±0 3 power-law den- 
sity profile (sometimes called the "central R~ 2 star clus- 
ter"; e.g., Becklin & Neugebauer 1968; Serabyn & Morris 
1995), and that a large-scale elliptical R~ 2 mass distri- 
bution around the GC is approximately consistent with 
the dynamical mass contained in the central 300 pc (see 
Haller et al. 1996). However, such a distribution cannot 
explain the observed distribution of light on scales larger 
than 30 pc, which is highly elliptical and cannot be de- 
scribed by a single power law. The observed NIR surface 
brightness distribution can much better be reproduced 
by a large-scale exponential disk with the Nuclear Stellar 
Cluster (NSC) at its centre. In the following, we call this 



disk the Nuclear Stellar Disk (NSD) and describe our two- 
component model. 

1. Nuclear Stellar Cluster: The radial K-band sur- 
face brightness distribution p h (R) in the central Rgc = 
0.625 pc (Paper I, Fig. 6, here corrected for extinction) can 
be best fitted by a volume density distribution of the form: 



p L (R) 



i± 



/ \ 



(io) 



with n = 2.0 and a core radius i?o = 0.22 pc. The ob- 
served K-band flux density and NIR luminosity due to 
cool stars in this area of A — 1 pc 2 are Sk ~ 350 Jy and 
icooi* ~ 3.1xl0 6 L Q , respectively (Paper I, Table 1). With 
these values we derive p% ~ 6.6 x 10 6 L Q pc~ 3 for Eq. ©. 
To fit both the observed large-scale surface brightness dis- 
tribution of the NB and the dynamical mass in the inner 
Rgc < 20 pc (see Sect. 5A and Fig. |lj), a steeper power- 
law dependence of the stellar density in this cluster is re- 
quired outside Rqc ~5-10pc. We obtained a good fit with 
n ~ 3 beyond 6 pc, but could not put strong constraints on 
where and how this transition in the slope exactly occurs. 
With these parameters we derive a total NIR luminosity 
of the NSC of Lnir = 6 ± 3 x 10 7 L . The relatively large 
uncertainty reflects the range of possible cut-off radii and 
outer slopes. The luminosity contribution from high-mass 
MS stars and the mass distribution in the NSC will be 
discussed in Sects. and Parameters of the NSC 
are summarized in Table 0. 

2. Nuclear Stellar Disk: Although other flat configura- 
tions may be conceivable for the NSD (e.g., a symmetric 
bar), we adopt a disk as the most likely configuration. The 
COBE NIR data are inconsistent with the disk having a 
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Fig. 12. Longitude and latitude profiles of the stellar NIR 
emission from the NB. Longitude and latitude offsets are 
measured from the position of SgrA*. a) and b) NIR 
volume emissivity profiles of the NB model representing 
the distribution of cool stars, c) and d) Integrated 4.9 /im 
surface brightness distribution of the volume emissivity 
model. For comparison, the thick grey line shows the IRAS 
60 fim latitude profile (90" resolution), scaled to match the 
4.9 /xm brightness, e) and f) Observed 4.9 /im (COBE) 
surface brightness profiles, corrected for extinction by dust 
in front of and inside the NB. Data points from positive 
and negative longitude and latitude offsets are plotted as 
filled and open circles, respectively. The average profiles 
were used to fit the models. The scatter of the data points 
indicates the longitudinal asymmetry, which we consider 
as the uncertainty of the model fit. The dashed profiles 
show the COBE beam. Solid curves show the profiles of 
the beam-convolved model map. 

single power-law emissivity profile. Instead, we obtain a 
good fit with an exponential-to-a-power law radial depen- 
dence of the volume emissivity: 



p h (R) oc exp 



mi 

2 



R 



a/2 



(11) 



with a the FWHM radius. Although the sum of two com- 
ponents with ax/2 — 120 pc and a 2 /2 = 220 pc, and 
n = 5 ± 1 fits the data best, the presence of two distinct 
components may not have a physical meaning, but could 
reflect extinction effects insid e the NB not accounted for 
in our model (see Sect. |5.5|) . The radial dependence of 
the volume emissivity in this disk can also be approxi- 
mated by a range of different power laws: flat (oc R~ 01 ) 
at R < 120 pc, oc R~ 3 - 5 for 120 pc < R < 220pc, and a 
steep drop-off (oc R~ w ) at R > 220 pc. Since the NB is 
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Fig. 13. Signatures of the distribution of massive stars 
in the Nuclear Bulge, a) Distribution of warm dust (50 
100 K, 60 Mm, IRAS), b) Distribution of hot dust (200- 
500 K, 12 /im, IRAS), c) Distribution of free-free (and 
non-thermal synchrotron) emission (6 cm, Parkes). The 
plot scales are logarithmic and cut-off levels are set to sup- 
press extended low-level emission. The most prominent 
Galactic centre features are labeled. Overlayed are con- 
tours of the NIR surface brightness of the Stellar Nuclear 
disk model (see Fig. O) . 



unresolved in latitude, only an upper limit for its FWHM 
scale height of 60 pc (0?40) could be derived from the 
COBE data. In Sect. |5.5| we derive a scale height of warm 
dust, which is closely related to the stellar NB, of 45±5pc. 
Adopting this scale height and n — 1.4±0.2 (see Sect. 5.5), 
Eq. ( |i~i"| ) yields a good fit to the latitude profile of COBE 
NIR data (Fig. |Jjf). The total NIR luminosity and aver- 
age midplane volume emissivity of the NSD then amount 
to Lnir = 6.5 ± 2 x 10 8 L© and ~ 100 L Q pc~ 3 , respec- 
tively. Parameters of the NSD are summarized in Table 



The corresponding total NIR surface brightness distri- 
bution, i.e., the sum of the contributions from NSC and 
NSD, is shown in Figs . |Tl"| b and[T^c and d. For comparison, 
we also show in Fig. |12| d the observed latitude profile of 
the IRAS 60 ^m emission from the NB (90" resolution). 
Finally, Figs. p"l| a and |l2| e and f compare, after convo- 
lution with the COBE beam, the modeled and observed 
NIR surface brightness profiles. Since we assumed symme- 
try with respect to the GC, the model was optimized to 
fit the average surface brightness at positive and negative 
latitude and longitude offsets, respectively. The observed 
and modeled latitude profiles do not match perfectly (Fig. 
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|l2|f). However, the small scale height of the NB together 
with the uncertainties in the beam shape, which may have 
slightly varied during the observations, do not allow to put 
better constraints on the vertical density profile of the NB. 

Although we cannot directly observe the emission from 
massive hot MS stars in the NB, their distribution can 
be derived from the emission by warm and hot dust and 
from free-free emission from HII regions. Figure [l3| com- 
pares the NIR morphology of the NB (cool stars, NIR 
model) with the distribution of the above signatures of 
massive stars. It becomes clear that massive stars are nei- 
ther evenly distributed throughout the NB nor are they 
present in the centre only. They are rather concentrated in 
several individual clusters, the brightest of which is cen- 
tered at Sgr A* and represents the centre of the the NSC. 
This confirms our earlier finding that young massive stars 
are strongly concentrated towards the centre (Paper I), 
but also shows that these stars are not exclusively formed 
in the central few pc of the NB. Our models account for 
this by extrapolating the (smooth ) ra dial distribution of 
the NSC to large radii (see Sect. 5.3). A comprehensive 
overview over the distribution of HII regions and non- 
thermal features in the GC Region can be found in LaRosa 
et al. (2000). 

5.3. The stellar luminosity of the Nuclear Bulge 

The total stellar luminosity of the NB depends on the con- 
tributions by high-mass stars whose optical/UV emission 
cannot directly be observed due to the high extinction. 
It has been shown in Paper II that this contribution de- 
pends critically on the Present-Day Mass Function which 
in turn is determined by the star formation history of the 
NB (continuous star formation, star bursts, etc.). In view 
of these uncertainties we use the integrated NIR SED as 
the only observable quantity to derive the total luminosity 
of cool stars in the NB and discuss possible corrections to 
this value due to contributions by high-mass stars. 

The total luminosity of cool stars in the NB can be de- 
rived by additionally correcting the observed stellar NIR 
emission (corrected for foreground extinction) for extinc- 
tion by dust inside the NB. Using the dust distribution 
and extinction model described in Sect. |3.5.2| , we derive 
Lnir(NB) - L coo i*(NB) - 6.8 x 1O 8 L (Table §, line 7). 

Assuming isotropic extinction and emission from the 
NB and assuming further that stellar emission is the only 
dust-heating source and that all stellar emission short- 
ward of A = 1 /im is completely absorbed by dust and re- 
radiated as thermal continuum emission, the total stellar 
luminosity of the NB (i.e., L CO oi* + -^hot*) should be equal 
to the observed total luminosity (i.e., LNiR,obs + ^dust) 
when corrected for all foreground extinction. This value 
amounts to L*(NB) = 6.7 x 10 8 L Q (Table §, line 6, last 
column), very similar to the intrinsic luminosity of cool 
stars derived above. However, to conclude that hot mas- 
sive MS stars do not contribute to the bolometric lumi- 
nosity of the NB is apparently wrong. Due to the highly 



inhomogeneous distribution of ISM in the NB, particularly 
the extreme dumpiness and very small scale height (see 
Sect. 5.5), only a small fraction of the short-wavelength 
stellar emission is actually absorbed and re-radiated (i.e., 
^dust •C ihot*) and the total luminosity may be consider- 
ably underestimated. 

Therefore, we use the intrinsic NIR luminosities of the 



NSC and NSD derived from our models (Sect. 5.2), as- 
sume that they represent the total luminosity due to cool 
stars, and try to supplement these values for the contribu- 
tion by hot massive stars. In Paper I (Table lb) we derive 
for the central 1.25 pc a luminosity ratio of hot (i.e., young 
massive) stars to cool stars (i.e. low-mass MS stars plus 
Giants and Supergiants) of ~30. According to our analy- 
sis of the KLF of the central 30 pc in Paper II, the ratio 
ihot*/icooi* drops to 4-5 when averaged over 30 pc. This 
behaviour can be approximately reproduced by assuming 
that £hot*/£ C ooi* = 30 holds for the entire NSC and that 
^nir = L to t holds for the entire NSD, i.e., no massive 
MS stars are associated with our model NSD. Note that 
NSC and NSD are two idealized model components which 
overlap spatially. The extrapolation of the NSC out to 
R = 230 pc actually provides L, ~ 5 x 10 8 L Q due to 
massive MS stars outside R > 30 pc, the region which is 



dynamically dominated by the NSD (see Sect. 5.4). This is 
approximately consistent with the number of massive star 
clusters and giant HII regions observed in the NSD (Fig. 
|l3| ) if each of these regions emits 5 — 10 x 10 7 L Q (e.g., 
Figer et al. 1999). The NIR luminosities of NSC and NSD 
(Sect. |J and Table @) convert then into a total stellar 
luminosity of the NB of L*(NB) = 2.5 ± 1 x 10 9 L Q . The 
large uncertainty is mainly due to the unknown contribu- 
tion by massive stars outside the central 30 pc. We adopt 
this value as the currently best representation of the total 
stellar luminosity of the NB. It is somewhat larger than 
the earlier, independent estimate of L(NB) ~ 1.4 x 10 9 L© 
by Cox & Laureijs (1988; see MDZ96, Table 5), which, 
however, did not account for the contribution by hot mas- 
sive stars. The stellar luminosities for different regions of 
the NB as derived from this model are summarized in 
Table || Table || summarizes the notations used for differ- 
ent physical components and spatial regions and their in- 
terrelation. Note that different physical components like, 
e.g., NSC and NSD, do spatially overlap, and that the 
characteristics of different spatial regions listed in Table |^ 
represent the corresponding average of the different model 
components in that region. 



5.4. The stellar mass of the Nuclear Bulge 

In the following we use M*/Lnir ratios observed in dif- 
ferent regions to determine M*/Lnir = f{R), then de- 
rive individual M*/Lnir ratios for NSC and NSD, and 
use these ratios to convert the NIR luminosity distribu- 
tions of NSC and NSD into a photometric mass distribu- 
tion of the NB. Note that only average M*/Lnir ratios 
are known for certain areas and that this ratio depends 
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Fig. 14. Enclosed mass in spheres of radius Rgc for the 
inner 500 pc of the Galaxy. Thin solid lines show the mass 
distribution of NSC and NSD as derived from our pho- 



tometric models (Sect. 5.2). The dotted line is the photo- 
metric mass distribution of the GB. The dashed horizontal 
line is for the central Black Hole (Eckart & Genzel 1998). 
The rightmost solid line (Mism) is for the total intestel- 
lar mass in the NB and the outer CMZ (Sect. pD| , The 
thick grey line represents the total enclosed mass, i.e., the 
sum of all these components. Data points with error bars 
refer to the following dynamical mass estimates: filled tri- 
angles: McGinn et al. (1989); open circle: Rieke & Rieke 
(1988); filled squares: Lindqvist et al. (1992); open tri- 
angle: Liszt & Burton (1978) and Burton & Liszt (1978) 
(adapted from Haller et al. 1996). The dashed curve (a) 
shows, for comparison, the power-law approximation to 
the enclosed mass by Sanders & Lowinger (1972) (Eq. 
(2.1) in MDZ96). The dotted curve (b) shows the pho- 
tometric mass distribution of the NB (without GB) ac- 
cording to our model when integrated through the entire 
NB in cylinders of radius i?cc (i- e -, the "observed" mass 
in a given area). The open square shows the total "ob- 
served" stellar mass of the central R = 15 pc derived from 
the KLF of the K-band mosaic (Paper II). 



critically on the ratio of mass-dominating low-mass MS 
stars to luminosity-dominating evolved star s, which is not 
constant throughout the NB (see Sect. 5.1). Here, we as- 
sume that NSC and NSD have different, but intrinsically 
constant M* / Lnir ratios and calibrate their values to re- 
produce the effective observed M*/Lnir ratios in different 
areas where both components overlap. 

The total mass in the central 1.25 pc or 30", 3.3 x 
10 6 M Q , according to Genzel et al. (1997), is dominated by 



the central Black Hole (A/ B h ~ 2.6 x 10 6 M ) , which leaves 
a total stellar mass of Af* ~ 7x 10 5 M Q in the central pc 3 . 



Assuming that the stellar mass distribution within the 
NSC has the same radial dependence as the NIR luminos- 
ity, the central mass density in Eq. ( |l0| ) becomes Pq 1 ~ 

3.3 x 10 6 M pc" 3 . This yields M*/£ N m ~ 0.5 which we 
assume to hold for the entire NSC. The total stellar mass 
of the NSC becomes then M*(NSC) = 3 ± 1.5 x 10 7 M Q . 

In Paper I we derive for the central 30 pc M*/Lnir ~ 1- 
The increase of M*/Lnir when averaged over larger ar- 
eas is consistent with the notion that the density of low- 
luminosity giants and low and intermediate-mass MS stars 
(spectral type 09 or later), which dominate the total stel- 
lar mass, decreases more slowly with increasing radius 
than that of K-luminous giants, supergiants, and early 
O stars (Paper I, Fig. 6). An effective M»/Lnir ratio 
of 1 for the central 30 pc can be reproduced by adopt- 
ing M*/£nir = 2 for the NSD , the same value as de- 
rived for the GB (see Sect. 4.2). This is consistent with 
the finding that the KLFs for low and intermediate-mass 
MS stars in Nuclear and Galactic Bulge are very similar 
(Paper II). Individual high-mass star- forming regions in 
the outer NB, which supposedly have a lower M»/Lnir 
ratio, are accounted for by extrapolating the NSC out 
to R = 230 pc (see Sect. 5J3). Assuming a 30% uncer- 
tainty for the M*/Lnir ratio, the NIR luminosity of the 
NSD (Table @) converts into a photometric stellar mass of 

1.4 ± 0.6 x 10 9 M Q . The masses of the two model compo- 
nents NSC and NSD and of different spatial regions of the 
NB are given in Tables [l] and ^, respectively. 

Figure |l4| shows the mass distribution in spheres of ra- 
dius i?GC around the GC derived from our volume emissiv- 
ity model of the NB (NSC + NSD) together with the above 
discussed M»/Lnir ratios. Also shown are the mass of the 
central Black Hole (Eckart & Gen zel 1998), the mass dis- 
tribution in the GB (sec Sect. 4.2), and the d istr ibution of 
ISM in the NB and the outer CMZ (see Sect. |1| here cor- 
rected for He and metals) . The mass distribution is domi- 
nated by the central Black Hole at Rgc < 2 pc, the NSC at 
2pc< i?cc < 30 pc, the NSD at 30pc< R GC < 400 pc, 
and the GB at -Rgc > 400 pc. This photometrically de- 
rived mass distribution in the central 500 pc with the 
Black hole included is compared to different dynamical 
mass estimates (for references see caption Fig. [li]). Both 
mass distributions agree very well. The presence of the 
NSD, which was inferred from the large-scale NIR surface 
brightness distribution of the NB, provides a logical expla- 
nation for the "kink" in the dynamical mass distribution 
at i?GC ~ 20 pc and the steep rise between Rqq ~ 20 pc 
and 100 pc. The power-law approximation by Sanders & 
Lowinger (1972) (Fig. [l| curve a; cf. Eq. (2.1) in MDZ96) 
also gives a reasonable estimate of the total enclosed mass 
in the inner 500 pc. However, it is inconsistent with the 
observed highly elliptical NIR surface brightness distribu- 
tion of the NB and, in detail, does not fit the dynamical 
mass distribution well. Note that these authors assumed 
a constant M* / -Lnir ratio for the entire GC region. Also 
shown in Fig. [14| is the stellar mass distribution contained 
in corresponding circular areas of radius i?cc when inte- 
grated throughout the entire NB (curve b). It reproduces 
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correctly the integrated stellar mass in the central 30 pc 
derived from an analysis of the KLF of the K-band mosaic 
(M*(30pc) - 1.2 x 10 s Mq, Paper II). Curve (b) cannot 
be directly compared to dynamical masses which refer to 
spherical volumes of radius i?cc- 

With M*(NB) = M*(NSC)+M*(NSD) = (1.4±0.6) x 
1O 9 M , the NB accounts for -1-2% of the total stel- 
lar mass in the GD and corresponds to — 10% of the 
mass of the GB (see Sect. [I7J). The average mass den- 
sity of stars in the NB (outside the central — 20 pc) 
is then p^ 1 ~ (200±50) M pc" 3 , which is more than 40 
times higher than the expected central mass density in 
the GB when its exponential volume emissivity/density 
profile (Eq. (H)) is interpolated into the centre. 



5.5. Morphology of Interstellar Matter in the Galactic 
Centre Region and its relation to the Nuclear 
Bulge 



In Sect. 4.3.3 we derived a total hydrogen mass of 6.8 x 
10 7 M Q for the central — 1 kpc and confirmed the well- 
known asymmetry in the mass distribution of ISM in the 
CMZ. A comparison of the distribution of stars (Figs. [Ti"| a 
and |l~5|c) and that of ISM (Figs. p~5| a and |l6|b) indicates 
that not all of the ISM in the CMZ is associated with the 
stellar NB. There are several indicators that allow to dis- 
tinguish between dust inside and outside the stellar NB: 

1. The stellar radiation density inside the NB is consid- 
erably higher than outside. Therefore, dust located inside 
the NB must have a higher temperature and specific lumi- 
nosity than dust outside the NB. Figures fL5| b and d show 
that both quantities are strongly enhanced in the central 
\l\ < 0?6 (~ 25 K and 3OL /M Q , respectively), indicat- 
ing the presence of a warm and luminous inner region 
(R < 100 pc) in the NB. At -0?6 > I > -1?9, both quan- 
tities decrease to — 21 K and 13-15 L q /Mq, respectively, 
indicating the presence of colder dust in the outer region 
of the NB. At I < — 1?9, the approximate outer "edge" of 
the NB, both quantities drop sharply and approach 16- 
17 K and 2-3L /Mq, respectively, values typical for dust 
located well outside the stellar NB, i.e., in the GB. 

2. At negative I, the longitude profiles of the total dust 
luminosity and the 4.9 /im surface brightness agree well 
(Fig. p~5[c) , indicating that the main dust-heating sources 
have the same spatial distribution as evolved stars, which 
dominate the NIR luminosity. Dust-heating mechanisms 
will be discussed in Sect. [|. In contrast, there is a relative 
dust luminosity excess at positive I. 

3. At I > 0?3 both the dust temperature and the specific 
dust luminosity drop sharply from values typical for the 
inner NB to values typical for ISM outside the NB (Figs. 
|l5| b and d). Together with the large peak at I ~ 0?9 of 
the hydrogen column density (see Fig. |l5|a) , this indicates 
the presence of large amounts of cold dust along the line 
of sight at positive longitudes. Hence we can attribute the 
dust luminosity excess at positive I to cold dust located 
outside the stellar NB. By subtracting this excess from the 




Fig. 15. a) Longitude profile of hydrogen column den- 
sity in the CMZ (see Fig. [Ttfb). b) Mass-averaged dust 
temperature in the CMZ. c) Total dust luminosity per 
pc 2 (solid line). Grey dots show the (dereddened) 4.9 /im 
surface brightness which represents the stellar luminosity 
(scaled to match L<j)- Dotted line at positive I: assumed 
intrinsic dust luminosity of the NB (see Sect. 5.5). d) 
Specific dust luminosity (Ld/Mn, solid line). Dotted line: 
assumed intrinsic specific dust luminosity of the NB (see 
Sect. [5~J). e) IRAS 60 and 100 /im surface brightness pro- 
files, representing warm dust. 



total dust luminosity map, we obtain the luminosity map 
of dust associated with the stellar NB (longitude profile 
shown in Fig. |l5|c). 



Since we showed in Sect. 5.2 that the distribution of 
stars (i.e., the density of the stellar NIR radiation field) 
is approximately symmetric with respect to the GC, we 
make the reasonable assumption that the specific lumi- 
nosity (L(i/Mh) of dust inside the NB is the same at pos- 
itive / as it is at negative I, where no cold dust intersects 
the line of sight, and 'symmetrize' the map of the specific 
dust luminosity. The mass distribution of warm luminous 
dust inside the stellar NB is then derived by dividing the 
corrected dust luminosity map (profile shown in Fig. [l5|c, 
dotted curve) by the corrected map of the specific dust 
luminosity (profile shown as dotted curve in Fig. EF 
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Low-luminosity dust outside the NB was filtered out from 
this map by setting a cut-off level for the specific dust 
luminosity at 5 L /M@ . The resulting hydrogen column 
density map of the NB is shown in Fig. |l6|d. This map is 
almost symmetric with respect to the GC. 

However, because the column density map was calcu- 
lated with the mass-averaged dust temperature, it is only 
a first-order estimate to the actual distribution of ISM in 
the NB. Figure p~5| b shows that the dust temperature in 
the inner part of the NB is considerably higher than in the 
outer parts. While we do not expect line-of-sight averaging 
effects at \l\ > 1°, the value of ~ 25 K at \l\ < 0?4 rep- 
resents certainly an average of cold dust in the outer NB 
and warm dust in the inner NB. To correct the iVn map of 
the NB (Fig. |l6|d) for this dust temperature gradient and 
to derive the three-dimensional density distribution, we 
transformed it back to a 240 /jm surface brightness map 
using an effective dust temperature map (profile in Fig. 
|l5|b) symmetrized in the same way as the specific dust lu- 
minosity (Fig. |l5|d). We then constructed and optimized 
a volume density model of ISM in the NB to fit, after 
convolution with the COBE beam, both the 'observed' 
240 /mi surface brightness and line-of-sight integrated ef- 
fective dust temperature. Although the sparsely resolved 
COBE maps do not allow to put strong constraints on 
the shape of the radial density profile, the flatness of the 
longitude FIR emission profile is not consistent with a ho- 
mogeneously filled, optically thin disk. It rather suggests a 
torus-like configuration. Since molecular line data indicate 
the presence of a "180-pc Molecular Ring" (see Sect. ^j]for 
discussion and references), the input configuration for our 
model was a 180-pc torus with Td = 21 K (see Fig. |l5|b) 
and an inner disk with a central density depression and 
Td > 25 K. The free parameters were the (radial) thick- 
ness of the torus, its density profile, and the radial density 
profile and temperature of the inner disk. The best fit to 
the data was obtained with: 

1. A torus with Ro = 180 pc, Td = 21 K, a radial density 
profile of the form given given by Eq. ( |TT| ) with n = 4, 
itfwhm = lOOpc (i.e., R m = 130pc and i? ut = 230pc), 
and 7ijj = 370 ± 50 cm" 3 , and 

2. an inner disk with R out = 120 pc, Td = (31 ± 3) K, and 
rig = 250±50ccm. The size and shape of the inner density 
depression could, of course, not be constrained. 

The presence of a sharp transition in the physical char- 
acteristics of ISM and/or a change of the interstellar ra- 
diation field at Rqc ~ 100 — 120 pc is also confirmed by 
the extremely sharp drop-off of the 60 and 100 /im emis- 
sion at \l\ ~ 0?7, which indicates a nearly complete ab- 
sence of warm dust (Td > 30 K) in the outer torus at 
R > 120 pc (see Figs. |l|e and |l7|c). The total hydrogen 
masses of the inner disk and outer torus are ~4x 10 6 M 
and ~ 1.6 x 1O 7 M , respectively. These numbers are com- 
piled in Table [|. The notations used for the different stellar 
and ISM components of the NB and their interrelation are 
summarized in Table |[ 

Since the NB is basically unresolved in latitude in the 
COBE maps, only an upper limit for its full scale height 



of 60 pc (0?40) could be derived. The IRAS maps provide 
a much higher angular resolution. The 100 /im emission is 
already representative for the bulk of the dust mass, al- 
though it may be somewhat biased towards warmer dust. 
A comparison between the COBE 60 /im, 100 fim, and 
240 /im maps shows that cold dust extends further than 
warm dust in longitude, but not in latitude. Since the dust 
ridge is very prominent in the 100 /im IRAS map and its 
scale height is considerably larger than the IRAS beam, 
saturation of the 100 /im IRAS map in the very centre does 
not falsify the result much. The vertical surface brightness 
profile could be best fitted by a volume density distribu- 
tion of the form of Eq. ( [Tl|) with n = 1.4 and a = 45 ±5 pc 

(0?30±0?03). The distribution is centered at b 0.05°, 

the position of SgrA* . This is in good agreement with the 
FWHM scale height of 12 CO at I = 0° of 50 pc (Bitran et 
al. 1997; see Fig. fl6|a). Since the 12 CO emission is optically 
thick, the actual scale height is expected to be somewhat 
overestimated by this tracer. 

In Fig. [l6|g, the resulting N^l map of the NB is dis- 
played, and Fig. [j~6]f shows the corresponding longitude 
profile as well as a comparison between the modeled 
(beam-convolved) and observed 240 /im surface brightness 
profiles. The corresponding radial profiles of the mid-plane 
hydrogen density and mass-averaged dust temperature are 
shown in Fig. |l7|. We call this configuration, i.e., the in- 
ner warm disk and the outer cold torus, the "Nuclear 
Molecular Disk" . It has the same size as the Nuclear 
Stellar Disk, i.e., i? ou t = 230 pc and A/ifwhm = 45±5pc. 

Note that the continuum data used in this paper do not 
allow to distinguish between a disk/ torus-like and a bar- 
like morphology. For example, the CMZ in IC342, which 
resembles in many aspects our own GC region, seems to 
resemble a bent bar rather then a disk/torus (e.g., Schulz 
et al. 2001). However, the main parameters of the NB 
(compiled in Table |^) would not change dramatically if 
a bar-like morphology is assumed rather than a disk or 
torus-like morphology. 

We attribute all other cold, low-specific-luminosity 
dust in the central kpc to ISM not heated by stars in 
the NB, and therefore located outside the stellar NB. 
Figure p~6| b shows the column density distribution of this 
material. The hydrogen masses at positive and negative 
I are ~ 2.9 x 1O 7 M and M£ ~ 1.1 x 1O 7 M , 

respectively. The average dust temperature and specific 
dust luminosity of this material are Td ~ 15 — 17 K 
and Ld/Mn ~ 2 — 3Lq/M q , respectively, compared to 
- 20 - 30 K and 10 - 5OL /M Q for dust inside the NB. 
The total hydrogen mass in the central 1 kpc of our Galaxy 
(the CMZ) amounts then to Mg MZ ~ 6 x 1O 7 M . Of 
this, one third is located in the Nuclear Molecular Disk, 
which is associated with the Nuclear Stellar Disk (i.e., 
the NB). Our estimate of the total hydrogen mass in the 
CMZ is intermediate between the values of 2 x 10 7 M 
and 2 x 10 s M derived by Mauersberger & Bronfman 
(1998) from the optically thin C 18 O(l-0) and optically 
thick 12 CO(2-l) lines, respectively, and agrees well with 
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Af H ~2-6x10 7 M q derived by Sodroski et al. (1995) from 
COBE DIRBE data. 

The total mass of the cold high-column density com- 
plex at I ~ 0?9 is Af H ~ 1O 7 M , of which the SgrB2 
GMC at I ~ 0?7 contains -10% (Gordon et al. 1993). 
If this complex would be located in front of the NB, it 
would cause high extinction. We assumed that the same 
low fraction of the total mass in this complex as inside 
the NB contributes to diffuse extinction while most of the 
mass is concentrated in small, ultra-opaque clouds which 
block 10% of the light (see Sect. |5.6| ), and tested the effect 
of this cloud complex on the final (extinction-corrected) 
stellar NIR maps of the NB. Figure [18] shows clearly that 
the model with the cloud complex located in front of the 
NB yields a strong, implausible 2.2 /mi excess, while the 
model with the cloud complex on the back side yields a 
good agreement of the dereddened 2.2 /im and 4.9 /mi sur- 
face brightness profiles. There is no reason to believe that 
the average effective stellar temperature in the NB is much 
higher at positive longitudes than at negative longitudes, 
nor does the specific dust luminosity indicate an excess of 
hot stars at I > 0.5° (see Fig. |l5|d). Since this test may 
be affected by the low angular resolution of the COBE 
data, we cannot ultimately conclude that the entire "+1° 
cloud complex" is located behind the NB. Indeed, the new 
2MASS and MSX images^ of the NB indicate the pres- 
ence of some very opaque dust lanes in front of the NB at 
I > 0.4°. However, on a large scale, we can exclude that 
this materi al causes strong extinction towards the NB (see 
Sect. |3.5.2| ). 



5.6. Physical characteristics of Interstellar Matter in 
the Nuclear Bulge 

From the model described in the previous section, we de- 
rive a total hydrogen mass of the Nuclear Molecular Disk 
of - 2 x 10 7 M Q , of which 4 x 10 6 M are located in the 
inner warm disk at R < 120 pc and 1.6 x 10 7 M Q (-80%) 
are distributed in the outer cold torus. The corresponding 
dust luminosities amount to — 2 x 10 8 L Q in each of the 
two components. Average densities in the inner disk and 
outer torus are (tih) ~ 100 cm -3 and ~ 150 cm" 3 , respec- 
tively. This high average hydrogen density supports the 
observation that most of the gas in the NB is in molecu- 
lar form. The mid-plane hydrogen column density through 
the entire NB is ~ 1 x 10 23 cm" 2 . Assuming homogeneous 
matter distribution, this corresponds to a total visual ex- 
tinction of ~ 100 mag between the front side of the NB 
and the GC, of which ~80mag are due to the cold outer 
dust torus. This clearly contradicts both the value of ~ 10- 
15 m ag deriv ed from the NIR colours of the NB (Sect. 
Lland 4.3.2| ) and the total visual extinction of 30 mag be- 
tween the Sun and the GC (e.g., Whitelock & Glass 1990). 



3 2MASS Image mosaic and 2MASS-MSX image combina- 
t ion by E. Kopan (IPAQ 



( |http:/ /antwrp.gsfc.nasa.gov/apod/ap000705.html[ ) 



This apparent contradiction can be solved when we 
assume that only ~10% of ISM (by mass) in the NB 
is homogeneously distributed, while 90% of the mass is 
trapped in small, discrete clouds which cover only a small 
fraction of the total area. A high degree of dumpiness 
is indeed confirmed by high-resolution submillimetre con- 
tinuum maps of the NB at A 800 /mi (Lis & Carlstrom 
1994) and A 1.2 mm (Zylka et al. 2001, in prep.). These 
maps show only very little extended emission between 
and around the dense clumps/clouds. We also find that 
the average column density of the central 12 pc o f A/h ~ 
7 x 10 22 cm" 2 (not deconvolved; see Sect. 4.3.3 ) agrees 
well with the value of ~ 8 x 10 22 cm" 2 derived from 
the A 1.2 mm map (when smoothed to the same angu- 
lar resolution and assuming Td = 25 K). It also agrees 
well with the value 7Vh ~ 8 x 10 22 cm" 2 which we de- 
rive from the C 18 O(l-0) map by Dahmen et al. (1997) 
using Eq. (7) of Mauersberger & Bronfman (1998) with 
[C 18 0]/[H 2 ] - 4 x 10" 7 and adopting T kin - 20 . . .30K. 
The fact that the chopped submm continuum maps, the 
total power low-resolution COBE FIR maps as well as the 
C 18 data see all the same hydrogen mass in the cen- 
tral 120 pc rules out that the high-resolution submm mea- 
surements have resolved out a considerable amount of ex- 
tended emission. This confirms the assumption that most 
of the ISM in the NB is concentrated in very compact 
molecular clouds. A first inspection of the 1.2 mm map 
suggests an area filling factor of dense clouds in the NB of 
~ 15% or less. This translates into a volume filling factor 
of only a few per cent. The average density of these com- 
pact molecular clouds, which contain ~ 90% of the total 
hydrogen mass in the NB, becomes then nn ~ 10 4 cm" 3 . 
Note that this is a very rough estimate and that the ac- 
tual densities in the cloud cores may be much higher. 
The average visual extinction of these clouds would be 
of order 1000 mag and they become optically thin only at 
A > 100 /mi. Thus, our mass estimate, which is based on 
the assumption that the 240 /mi emission is optically thin, 
should still be valid. We call these clouds "ultra- opaque" 
and assume that ~ 10% of the stellar emission from the 
NB and the back si de of the GB is completely blocked out 
by them (see Sect. |3.5.2 ). Further characteristics of these 
clouds are being investigated in Paper IV (Zylka et al., in 
prep.). 

The physical reason for this dumpiness is most prob- 
ably the tidal stability limit. We find that in the gravita- 
tional potential of the mass distribution shown in Fig. 
[il] only clouds with densities higher than a few times 
10 3 — 10 4 cm -3 are stabilized against tidal forces by their 
own gravitation (cf. Glisten & Downes 1980). In contrast 
to Giisten & Downes (1980), who use the mass distribu- 
tion derived by Sanders & Lowinger (1972), we find that 
this limit holds only for the outer NB (R > 120 pc) and 
that the shallow mass distribution in the inner Nuclear 
Stellar Disk causes a stability window between i?GC ~ 50 
and 110 pc. Inside 50 pc only superdense clouds are stable 
against tidal disruption. 
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For the central region of IC342, which resembles in 
many aspects our own GC region, but is seen face-on, 
Schulz et al. (2001) also hnd that the molecular gas is 
extremely clumpy and that the bulk of the gas mass in 
the inner — 50pc of this galaxy is concentrated in com- 
pact clouds with A/h ~ 1O 6 M and nn ~ 10 4 cm~ 3 . 
Quantitatively, this is a striking similarity to the condi- 
tions in our own GC region. 

The total mass and average density of the homoge- 
neously distributed intercloud medium in the NB are 
then M H = 2.0 ± 0.5 x 10 6 M Q and n H ~ lOcnT 3 , re- 
spectively. The presence of a diffuse, hot, ionized inter- 
cloud medium in the NB is supported by the presence 
of extended thermal (free-free) radio emission (see Fig ^, 
6cm map). Mezger & Pauls (1979) derived for the ex- 
tended, low-density ionized gas in the NB an average elec- 
tron temperature of T c ~ 5000 K and an electron density 
of n e ~ 10 cm -2 . The total mass of this component is 
Mhii ~ 2 x 10 6 Mq. Both values are in good agreement 
with our estimate. The detection of extended X-ray emis- 
sion in the inner ~ 1?8 of the NB indicates the presence 
of second diffuse component, i.e. a very hot Plasma (e.g., 
Koyama et al. 1996). 



6. Towards an empirical model of the Nuclear 
Bulge 

Using COBE and IRAS NIR to FIR data, we derived the 
morphology and physical characteristics of stars and ISM 
in the central kpc of the Galaxy. Figure |l9| summarizes 
the projected large-scale distribution of stars and ISM 
in the Galactic Centre Region. We find that the central 
D ~ 500 pc are dominated by the Nuclear Bulge, a dense 
and massive disk-like complex of stars and ISM, which is 
clearly distinguished from the much larger Galactic Bulge 
in which it is embedded. Similar distinct nuclear star clus- 
ters have been found by HST observations in the centers 
of many exponential bulges in intermediate and late-type 
spiral galaxies (e.g., Carollo et al. 1998). Here, we derived 
density and luminosity distributions of the various stellar 
and gaseous components populating the Nuclear Bulge of 
our own Galaxy. Notations and spatial extents of these 
components are described in Table pi and their physical 
characteristics are compiled in TablelS. 

The Nuclear Bulge appears, in projection, as a flat bar 
with an outer radius of 230±20pc, FWHM scale height 
45±5pc, and diameter-to-thickness ratio 5:1. Its overall 
morphology appears disk-like and symmetric with respect 
to the Galactic Centre. This disk may be elliptical, in 
which case it would rather resemble a bar. The Nuclear 
Bulge consists of the large Nuclear Stellar Disk with the 
R~ 2 Nuclear Stellar Cluster at its centre, and the Nuclear 
Molecular Disk. In our analysis, we neglected the small- 
scale asym metr y seen in the COBE NIR images (see Fig. 
and Sect. ^2). On small scales, the Nuclear Bulge may 
well be not as symmetric as suggested by our results. This 
remains, however, subject to further studies which have to 



combine NIR, MIR, and submm data with much higher 
angular resolution than provided by the COBE data. 

The total mass of the Nuclear Bulge amounts to 
1.4±0.6x 10 9 M Q , of which ~ 99% is contained in stars and 
1% in ISM, mainly in the form of dense molecular clouds. 
The Nuclear Bulge accounts thus for ~1% of the total 
stellar mass in the Galactic Disk and Bulge. In contrast 
to the usual assumption that the mass distribution in the 
central 300 pc is dominated by an R~ 2 cluster, we find that 
the density profile of the Nuclear Stellar Cluster steepens 
outside i?cc ~ 5 — 10 pc and that it only dominates the 
inner Rgc ~ 20 — 30 pc. At larger radii, the mass distribu- 
tion is dominated by the Nuclear Stellar Disk which has a 
rather flat radial density profile inside Rgc ~ 120 pc and 
a well-defined outer boundary at —230 pc. The Galactic 
Bulge begins to dominate the mass distribution only for 
Rgc > 400 pc. The photometric mass distribution in the 
inner 500 pc derived from our models of the Nuclear and 
Galactic Bulge (with the central Black Hole included) 
agrees well with the dynamical mass distribution derived 
by other authors and with the total stellar mass in the 
central 30 pc as deduced in Paper II from the KLF of this 
region (see Fig. |l4|). The total stellar luminosity of the 
Nuclear Bulge amounts to L* ~ 2.5 ± 1 x 10 9 Lq, of which 
7 ± 2 x 10 8 L Q (~30%) arise from cool low-mass MS and 
evolved stars. The contribution by hot and massive stars 
(~70%) is highly uncertain. The Nuclear Bulge thus con- 
tributes ~5% to the total luminosity of the Galactic Disk 
and Bulge. 

The stellar population of the Nuclear Bulge appears 
to resemble that of the Galactic Bulge. Unique for the 
Nuclear Bulge, however, is a population of young massive 
stars and an over-abundance of luminous giant and super- 
giant stars, which are both strongly concentrated towards 
the centre and are mainly associated with the Nuclear 
Stellar Cluster. However, the distribution of hot dust, HII 
regions, and supernova remnants (Fig. [l3|; see LaRosa et 
al. 2000) indicates that star formation occurs throughout 
the entire Nuclear Bulge and not exclusively at its centre. 
Ongoing star formation is one of the key characteristics 
which distinguish the Nuclear Bulge from the Galactic 
Bulge. The star formation history of the Nuclear Bulge, 
based on the analysis of the KLF of the central 30 pc, has 
been discussed in Paper II. 

The total hydrogen mass in the central kpc, referred 
to as the Central Molecular Zone, amounts to Mh ~ 
6 x 10 7 M or, if corrected for He and metals, to A/ism ~ 
1 x 10 8 M Q . Wc find that -2/3 of this material is lo- 
cated outside the stellar Nuclear Bulge and not heated 
by these stars. About 3 x 10 7 M Q are at positive I (0?5- 
4°) and 1 x 10 7 M Q at negative I (-2° --4°), thus caus- 
ing the observed asymmetry in the distribution of ISM in 
the Central Molecular Zone. The physical reason for this 
asymmetry is not known. Only — 2 x 1O 7 M are directly 
associated with the stellar Nuclear Bulge. Gas and dust 
in this Nuclear Molecular Disk, which is approximately 
symmetric with respect to the Galactic Centre, are dis- 
tributed in a thin and warm inner disk (Td — 30 K) and 
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Stellar Components: Explanation 

Nuclear Stellar Cluster (NSC): Central p oc r~ 2 stellar cluster, as first described by Becklin & Neugebauer (1968). Its 

central part is identical with the star cluster discussed in Paper I (Fig. 6). 

Nuclear Stellar Disk (NSD): Large-scale disk (R ou t ~ 230 pc) of mainly cool stars with a possible density drop at 

R ~ 120 pc. Contains most of the mass in the NB. 

Nuclear Bulge (NB): The sum of NSC and NSD. 

ISM Components: 

Inner warm disk: Warm dust and molecular gas in the NB which is restricted to R < 120 pc. 

Outer cold torus: Colder dust and molecular gas in the outer NB (still heated by these stars), possibly 

related to the "180 pc Molecular Ring". The reason for the sharp transition in ISM 

characteristics at R ~ 120 pc is not yet understood. 
Nuclear Molecular Disk (NMD): The sum of the inner warm disk and the outer cold torus. Represents all molecular gas 

which is physically associated with the stellar NB and heated by these stars. 
Central Molecular Zone (CMZ): All molecular gas in the central kpc of the Galaxy. Includes the Nuclear Molecular Disk 

and cold molecular gas outside the NB in the GB. 

Spatial regions: (see Table |^) 

Central 1.25 pc: Sphere of R = 0.625 pc around the GC (Volume ~ lpc 3 ). Dominated by the NSC. 

Inner NB: Stars and ISM at R < 120 pc, defined by the sharp transition of ISM characteristics 

and a possible drop of the stellar density at R ~ 120 pc. Includes the central 1.25 pc 
Outer NB: Outer part of the NB and NMD (the outer cold ISM torus), mainly defined by the 

complete lack of warm and hot dust at R > 120 pc (the outer cold ISM torus). 

Stellar population dominated by the NSD. 



a colder, dense outer torus (~20K) between Rqc ~ 120 
and 220 pc, which contains ~ 80% of the total gas mass 
in the Nuclear Bulge. 

The presence of such a molecular torus in a large-scale 
bar potential (due to the barred Galactic Bulge) is pre- 
dicted by the same dynamical models as described in Sect. 
[|. Due to shocks and angular momentum loss on the in- 
ner self-intersecting X\ orbits in the HI Central Disk, the 
gas is driven further inward and being compressed to den- 
sities at which molecular clouds form, and finally settles 
on stable X 2 orbits (e.g., Binney et al. 1991; Englmaier 
& Gerhard 1999; Englmaier & Shlosman 2000). The outer 
Central Molecular Zone would then represent the region 
between the innermost stable X\ orbit and the more cir- 
cular X2 orbits. The massive outer torus of the Nuclear 
Molecular Disk can be associated with the region where 
this gas has accumulated on stable X 2 orbits and becomes 
dense enough to form stars. One implication of this in- 
terpretation is that star formation in the Nuclear Bulge 
should occur - apart from the very centre - mainly in the 
outer molecular torus. Indeed, the distribution of giant 
H II region complexes and supernova remnants across the 
entire Nuclear Bulge, with the outermost ones being at 
Rgc = 170 - 190 pc (SgrD and E; see LaRosa et al. 2000 
and Fig. |l3|) is consistent with the presence of a modest 
star-forming ring with a mean radius of 180 pc around the 
galactic centre. The relation between this massive torus 
and the narrow 180-pc Molecular Ring (Bally et al. 1987) 
is not fully understood (see Sect. |l|). 



Extinction arguments lead to the conclusion that ISM 
in the Nuclear Molecular Disk is extremely clumpy. This 
is confirmed by high-resolution submm continuum maps 
of the inner Nuclear Bulge (Lis & Carlstrom 1994; Zylka 
et al., in prep.). We find that more than 90% of the total 
gas mass in the Nuclear Molecular Disk is concentrated in 
very compact, ultra-opaque molecular clouds with average 
densities of (nn) > 10 4 cm~ 3 and a volume filling factor of 
only a few per cent. These clouds are completely optically 
thick for UV, optical, and NIR radiation, but they become 
optically thin for the thermal FIR/submm radiation at 
wavelengths longer than 100 /im. The physical reason for 
this dumpiness is most probably the tidal stability limit, 
since only clouds with densities higher than a few times 
10 3 — 10 4 cm -3 are stabilized against tidal forces by their 
own gravitation in the gravitational potential of the mass 
distribution in the Nuclear Bulge. The remaining < 10% 
of interstellar mass in the Nuclear Bulge are contained in 
a diffuse, thin and warm intercloud medium. 

The dominant heating source of the thin intercloud 
medium in the Nuclear Bulge is the strong UV radiation 
field due to young massive stars in the Nuclear Stellar 
Cluster and other regions of high-mass star-formation in 
the outer Nuclear Bulge, as well as a widely distributed 
population of intermediate-mass MS stars. Due to the ex- 
treme dumpiness of the ISM, at least the soft UV radi- 
ation can penetrate the entire Nuclear Bulge, ionize the 
thin intercloud gas, and heat the dust in this compo- 
nent to high temperatures. This soft UV radiation to- 
gether with the strong stellar NIR radiation field also 
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Table 9. Adopted characteristics of the Nuclear Bulge 



Component: 


Total 


Central 1.25 pc 


Inner NB 


Outer NB 


Ref./Se 


Size: 














R/pc 


230±20 


0.625 


120±20 


120-230 


5^2 






/ifwhm/pc 


45±5 


- 


45±5 


45±5 


5.2 




5.5 


Stars: 














M*/M 


f .4 ± 0.6 x fO 9 


7 x 10 5 


8 ± 3 x 10 s 


6 ± 2 x 10 s 


5.4 




( P M t >/(M pc~ 3 ) 


200 ± 50 


7 x 10 5 


400 ± 100 


120 ± 30 


5.4 






7 ± 2 x f s 


1.4 x 10 6 


4 ± 1 x 10 s 


3 ± 1 x 10 s 


5.3 






2.5 ± f x f 9 


4 x 10 7 


2 ± 1 x 10 9 


5 ± 2 x 10 s 


5.3 




-k*,NIR/ -t^*,tot 


U.Zo 


U.UoO 


u.zu 


n f\n 


5.3 




(or * J/(Lm DC -3 ) 


fOO ± 50 


4 x 10 7 


200 ± 50 


60 ± 20 


5.3 




-t *,eff/ 


a Ann 4- 4nn 


>4500 


4400 4- 400 


44nn + 4nn 


4.3.2 




/ 71 iT IT \ /VA/T /T 

(M*/L*,nir)/(M©/L )^ ; 


2 ± 0.3 


0.5 ± 0.2 


2 ± 0.3 


2 ± 0.3 


5.3 




5.4 


(M,/L*,tot)/(M Q /L Q ) (i,) 


0.6 ± 0.3 


0.02 ±0.01 


0.4 ± 0.2 


1.2 ± 0.5 


5.3 




5.4 


ISM: 














Mh/M s 


2.0 ±0.3 x f0 7 


_(d) 


4 ± 1 x 10 6 


1.6 ±0.3 x 10 7 


5.5 




{ri H )/cm- 3(c) 


f 40 ± 30 


_(d) 


100 ± 30 


150 ± 30 


5_6 




idust/L 


4.0 ±0.2 x fO 8 




2 ± 0.3 x 10 s 


2 ±0.3 x 10 8 


5.5 




Tdust/K 


20-30 


_(d) 


31 ±3 


21 ± 1 


5_5 




(£d«t/M H )/(L Q /M Q ) 


20 ±3 


_(d) 


50 ± 10 


12.5 ± 2 


5.5 




ISM/Stars: 




_(d) 










(Mh/M*) 


1.4 ±0.5 x 10~ 2 


5 ±1.5 x 10 -3 


2.7 ±0.5 x 10~ 2 






(idust/^.NIR) 


0.6 ±0.2 


_(d) 


0.5 ±0.2 


0.7 ±0.2 






(£duSt/£*,tot) 


0.16 ±0.08 


_(d) 


0.10 ±0.05 


0.4 ±0.2 







(<0 



Tnir derived from NIR SEDs. Represents ~ T c g of cool stars and does not include contribution from hot massive stars 
Includes contribution from hot massive stars 

Average riH, assuming homogeneous matter distribution. The ISM is actually composed of discr ete dense molecular clouds 
with (nu) ~ 10 4 cm -3 and a thin intercloud medium with nu ~ 10 cm -3 (see discussion in Sect. |5.6[ ). 
Values could not be derived 

Note that some values for the three sub-regions of the NB represent average characteristics of different physical components 
which overlap in these regions. Some of these values are not individually derived in the corresponding sections, but arise 
from the combined model of the NB. 



dominates the heating of the dense molecular gas (e.g., 
Poglitsch et al. 1991; Hiittemeister et al. 1998; Schulz et 
al. 2001). The relative lack of hot dust (200 - 400 K) and 
the abrupt decrease of specific dust luminosity beyond 
i?GC ~ 100 — 120 pc supports our finding that most of the 
gas and dust in the Nuclear Molecular Disk is contained 
in very compact and massive molecular clouds distributed 
mainly in a torus between Rgc ~ 120 and 220 pc. Most 
of their mass then is effectively shielded from the strong 
interstellar radiation field in the Nuclear Bulge. Gas in 
the inner Nuclear Molecular Disk must also be clumpy so 
that the radiation field can penetrate and maintain the 
high mass-averaged dust temperature of > 30 K. However, 
the absence of a stringent tidal stability limit between 
-Rgc ~ 50 and 110 pc allows smaller clouds to exist in 
this region, which have a larger surface-to-mass ratio and 
can be more effectively heated, thus explaining the large 
excess of hot dust compared to the outer torus. 



With its flat disk and ring-like morphology, high stel- 
lar density, high degree of dumpiness of ISM, and ongoing 
star formation, the Nuclear Bulge may be tentatively con- 
sidered as counterpart of the central complex in IC 342 
(e.g., Wright et al. 1993; Boker et al. 1997; Mauersberger 
& Bronfman 1999; Schulz et al. 2001) or even of the mas- 
sive gaseous tori seen in starburst galaxies (e.g., Downes 
& Solomon 1998). 

Since we used a distance of 8.5 kpc to the Galactic 
Centre, all sizes have to be scaled by (8/8.5)^ 0.94, and 
masses and luminosities by (8/8. 5) 2 ~ 0.89, when the 
numbers are being compared to other studies which use a 
distance of 8 kpc. These corrections are, however, smaller 
than the uncertainties of our results. 
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Appendix A: The interstellar dust extinction curve 



Table A.l. Dust opacities (per gram of dust) 



A 

H 


[cm 2 g x ] 


A 

H 


k(A) 
Kg" 1 ] 


A 

H 


«(A) 
Kg" 1 ] 


0.55 


32100 


4.9 


926 


140 


20 


1.25 


9030 


12 


1190 


240 


7.3 


1.65 


5300 


25 


415 


800 


0.71 


2.2 


3590 


60 


77 


1300 


0.31 


3.5 


1800 


100 


34 


3000 


0.072 



A "standard" extinction curve was adopted for the 
entire GD and the NB since extinction along the line 
of sight toward the GC is believed to be dominated by 
dust in the diffuse ISM which accounts for the intercloud 
gas as well as for low- and intermediate-density molecu- 
lar clouds (Tielens et al. 1996). The dust opacity spec- 
trum used here is based on data from Rieke & Lebofsky 
(1985) for A = 0.36... 13 /xm and Draine & Lee (1984) 
for A = 13 /im. . . 1.3 mm. Tabulated values for dust grains 
without ice mantles and no coagulation were taken from 
Ossenkopf & Henning 1994, who assumed a M RN grain 
size distribution (Mathis et al. 1977). Figure A.l shows the 
dust opacity spectrum used in this paper and Table A.l 
lists opacities per gram of dust for selected wavelengths. 
A "standard" H-to-dust mass ratio in the solar neighbour- 
hood of (Mn/Md)o = 110 was adopted (Draine & Lee 
1984; Mathis & Wallenhorst 1981; Kriigel et al. 1990), 
which scales with the reciprocal of the relative metallicity 
Z/Zq (see Fig. B.l). The resulting Ay/Nn ratio amounts 
to (A Y /N H ) Q = 5.3 10~ 22 (for Z/Z Q = 1; Bohlin et al. 
1978). The FIR/submm spectral index of the dust opac- 
ity is (3 = -1.8. 



Appendix B: Distribution of Interstellar Matter in 
the Galactic Disk 

To treat the emission and extinction by dust in the GD 
in an appropriate way and to correctly derive NIR surface 
brightness maps of the NB, a three-dimensional model of 
the distribution of ISM in the GD was developed. The 
primary goal of this model was to derive extinction maps 
towards the GB and NB and to separate the MIR to FIR 
emission from the GD into contributions from the front 
and back sides with respect to the GC region. Proper 
modeling of the NIR surface brightness distribution of the 
GB, which is impossible without applying correct extinc- 
tion corrections, is particularly important since due to its 
exponential profile the separation of GB and NB is very 
sensitive to how the surface brightness of the GB is inter- 
polated into the region occupied by the NB. 

Here, the GD was modeled by an axisymmetric disk 
with a central hole and an overlayed spiral arm pattern. 
Three gas components mixed with dust were considered: 
i) atomic hydrogen (HI), ii) molecular Hydrogen (H 2 ), 
and Hi) HII. The dust model is described in Appendix 
[a|. Radial density profiles at b = 0° were adopted from 
Liszt (1992) for HI and from Bronfman (1992) for H 2 , 
and were adjusted to fit the FIR data. A sech 2 vertical 
density profile was adopted for the H 2 disk and a double- 
Gaussian profile for the H I disk. The scale heights of the 
H I and H 2 disks as function of galactocentric radius - in- 
cluding the exponential flare and the warp of the outer 
disk - were adapted to various measurements summarized 
by Wouterloot et al. (1990). The dust temperatures were, 
to first order, adopted from Sodroski et al. (1994) and then 
varied to fit the COBE FIR data with the above model. 
The metal abundance gradient was taken from Afflerbach 
et al. (1997). For the GC region we use a relative metal- 
licity of Z/Zq = 2 as derived by Mezger et al. (1979). 
A four-arm, modified logarithmic spiral density pattern 
was overlayed on the axisymmetric disk. The general spi- 
ral geometry was adopted from Georgelin & Georgelin 
(1976) and Englmaier & Gerhard (1999), and adjusted to 
fit the COBE FIR data. It resembles the pattern derived 
by Drimmel & Spergel (2001), but is not identical. The de- 
tails of this spiral arm model may not necessarily reflect 
the true spiral pattern of the Galaxy. However, it fits the 
COBE FIR data sufficiently well for our purpose and pro- 
vides a reasonable estimate of how the extinction towards 
the GB rises with increasing galactocentric radius towards 
the GD molecular ring. The spiral arms have a Gaussian 
density cross section. The density contrast between spiral 
arms and inter-arm regions and its radial gradient were 
adopted from Heyer & Terebey (1998). The local Orion 
arm was not included in our model. The main parameters 



of the GD model are summarized in Fig. B.l 

With these parameters, the GD was modeled ou t to 
R = 20kpc with a grid cell size of (100 pc) 3 . Figure B^ 
compares the observed 240 /im emission profiles with the 
profiles predicted by the model (see Fig. || for the latitude 
extinction profile). The GD model contains Mh 2 ~ 1 X 



26 



Launhardt et al.: Large-scale characteristics of the Nuclear Bulge 



10 9 M Q molecular hydrogen and Mm ~ 3 x 10 9 M Q atomic 
hydrogen of which ~ 2 x 10 9 M Q are located at R < 15 kpc. 
These values are in good agreement with other estimates. 

This model allows to derive the thermal emission and 
extinction due to dust along any line of sight in the inner 
Galaxy. Although our model is not intended to represent 
the exact morphology of the entire GD and does not in- 
clude the stellar disk, its strength is that it is based on 
a complete set of parameters derived from observations 
without artificial assumptions and that it reproduces the 
observed FIR emission and colour temperatures and ex- 
tinction in the inner Galaxy. A much more sophisticated 
three-dimensional model of the GD. which is also based 
on COBE/DIRBE data, has recently been developed by 
Drimmel & Spergel (2001). 
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Fig. 16. a) Distribution of 12 CO(1-0) (HPBW = 0?15; 
Bitran et al. 1997; see also Mauersberger & Bronfman 
1998). b) Hydrogen column density map of the cen- 
tral 8° derived from the COBE 240 /im dust emission 
map (HPBW = 0?7). Contributions from the GD are sub- 
tracted. Contour lines are at 1.5, 2.5, 3.5 to 15.5 by 
2 x 10 22 cm~ 2 . c) Column density map of ISM outside 
the NB. d) Column density map of ISM in the NB. e) 
Hydrogen column density map of the NB model (see Sect. 
|5.5| ). Contours are in steps of 10% of the maximum. Dotted 
contours show the distribution of ISM outside the NB. f ) 
Longitude profile of hydrogen column density inside (black 
solid line; cf. map e) and outside the NB (dotted curve, 
not deconvolved; cf. map c). Light grey curve: observed 
240 fim surface brightness profile of dust inside the NB. 
Dashed line: modeled 240 /im surface brightness profile of 
dust inside the NB (convolved with the COBE beam). 
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Fig. 17. a) Modeled radial mid-plane profile of the aver- 
age hydrogen number density in the NB (assuming ho- 
mogeneous matter distribution; see discussion in Sect. 
5.5). Note that the central density depression is not well- 
constrained, b) Mass-averaged dust temperature, c) 
IRAS 60 and 100 /im surface brightness profiles, represent- 
ing warm dust (positive and negative offsets from the GC 
averaged) . 
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Fig. 18. Longitude profiles of the 2.2 (im and 4.9 /im 

emission from the NB (COBE). The observed emission 
is dereddened with two different extinction models, as- 
suming the "+1° cloud complex" being located in front 
(dashed curve and empty circles) and behind the stellar 
NB (solid curve and filled circles). 
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Fig. 19. a) Model of the completely dereddened 4.9 /im 
surface brightness of the NB. b) IRAS 60 /im emission 
from warm dust, c) Model of the hydrogen column den- 
sity in the NB, representing mainly the distribution of 
cold dust. Dashed lines: Ah contours of cold ISM outside 
the NB (as observed with the large COBE beam). d) 
Schematic overview of the Galactic Centre Region viewed 
from the Sun. Dotted lines: stellar NIR emission from 
the bar-like Galactic Bulge. Solid lines: cold ISM in the 
Central Molecular Zone, located outside the NB. Grey- 
scale plot: the Nuclear Bulge. 
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Fig. A.l. Opacity spectrum of the dust model used in 
this paper (cm 2 per gram of dust). 
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Fig. B.l. Main parameters of the Galactic Disk model, 
a) HWHM height of the density profiles of H I and H2 in 
the GD. Data points (black squares for H2 and grey circles 
for HI) were taken from Wouterloot et al. (1990). The 
solid lines show the scale heights of H I and H2 adopted 
for the model, b) Ring-averaged volume density profiles 
( H cm" 3 ) at b = 0° The ripples are due to the spiral 
structure, c) Mean dust temperature profiles of the HI 
and H 2 components of the Galactic Disk. d) Relative 
abundance of heavy elements as function of galactic radius 
(after Afflerbach et al. 1997). 
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Fig. B.2. Observed and modeled 240 /im emission pro- 
files. Grey squares and thick grey curve: COBE data, 
Black solid lines: model, a) Average latitude profile of 
the Galactic Disk between — 20° < I < 20°, excluding spiral 
arms and the NB region, b) Longitude profile at b = 0°. 
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